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Abstract 

We review the prospects for bottom production physics at the LHC. 

^ ; 1. INTRODUCTION 

In the context of the LHC experiments, the physics of bottom flavoured hadrons enters in different 
^ ■ contexts. It can be used for QCD tests, it affects the possibilities of B decays studies, and it is an 

Tj- . important source of background for several processes of interest. 



The physics of b production at hadron colliders has a rather long story, dating back to its first 



observation in the UA1 experiment. Subsequently, b production has been studied at the Tevatron. Besides 
the transverse momentum spectrum of a single b, it has also become possible, in recent time, to study 
correlations in the production characteristics of the b and the b. 

At the LHC new opportunities will be offered by the high statistics and the high energy reach. One 
expects to be able to study the transverse momentum spectrum at higher transverse momenta, and also 
[ to exploit the large statistics to perform more accurate studies of correlations. 

This chapter is organized as follows. 

Section ^| is mostly theoretical. Its goal is to provide benchmark cross sections and distributions 
for the LHC, including rates relevant for the trigger requirements of the experiments. Furthermore, a 
discussion of the present status of b production phenomenology at hadron colliders is given. In this 
context, one cannot forget that the theoretical status is a mixed success. On one side, the shape of 
distributions and correlations are reasonably well explained by perturbative QCD. On the other side, 
however, the observed cross sections at the Tevatron are larger than QCD predictions. It is hoped that 
further studies may help to understand the nature of the discrepancy. As of now, we see two possible 
explanations: either the absolute normalization of the cross section is not correctly predicted due to the 
presence of large higher order terms, or the shape of the distributions is distorted by some perturbative 
or non-perturbative effects (like, for example, fragmentation effects). With the wide p T range covered 
by the LHC experiments, and perhaps also with the possibility of performing more accurate studies of 
correlations, these two possibilities may be distinguished. The problem of fragmentation effects has been 
studied in this workshop also from the point of view of hadronization models in Monte Carlo programs, 
in Section ||. This study deals with the hadronization model in the HERWIG Monte Carlo program. 
Its aim was to understand whether, in simple realistic models of hadronization, the usual assumption 
of QCD factorization is really at work. In general, the problem of studying how realistic is the heavy 
flavour production mechanism implemented in shower Monte Carlo's is quite important, and probably 
will require a considerable effort. Along this line, in Section |], a problem in the heavy flavour production 
mechanisms implemented in Pythia is examined. 

Further self-contained theoretical topics are dealt with in Sections || and ||. Section || deals with 
the charge asymmetry in b production in pp collisions. In this context, QCD is not of great help, since in 



perturbative QCD charge asymmetries turn out to be extremely small. Instead, studies are made within 
specific hadronization models, that are parametrized in such a way that they fit charm asymmetry data. 
This topic, besides being interesting in its own, since it deals with a phenomenon which is dominated by 
non-perturbative physics, has also an impact on CP violation studies in B decays. 

In Section g quarkonium production is discussed. This subject has been intensively studied in 
recent years, following an initial CDF observation of a J/* production rate much higher than theoretical 
predictions. This has triggered, from the theoretical side, the understanding that the fragmentation pro- 
cess is the dominant mechanism in quarkonium production. Besides this, a novel branch of applications 
of perturbative QCD, the NRQCD approach, has emerged, that may be useful to explain the production 
process. 

In Section [7|, the prospects for b detection are discussed. It is shown that there is a complemen- 
tarity between ATLAS/CMS and the LHCb experiment, with a certain region of overlap. In particular, 
the LHCb experiment can detect very low momentum heavy quarks, while the other experiments can 
reach the very high transverse momentum region. Some results on correlations measurements are also 
given, exploring the possibility of looking at one b decaying into a J/^, and the other decaying semilep- 
tonically. Double heavy flavour production, charge asymmetry, polarization effects, and doubly-heavy 
meson production are also discussed. 

In Section |8| the tuning of the multiple interaction parameters in Pythia is illustrated. The correct 
treatment of multiple interactions is important to model the multiplicity observables in both minimum- 
bias and heavy flavour events. 

2. BENCHMARK CROSS SECTIONS^ 
2.1 Total cross sections 

It is assumed that heavy flavour production in hadronic collisions can be described in the usual improved 
parton model approach, where light partons in the incoming hadrons collide and produce a heavy quark- 
antiquark pair via elementary strong interaction vertices, like, for example, in the diagram of fig. [I]. This 




Fig. 1: Typical diagram for heavy flavour production 



description is appropriate for all hard processes in hadronic collisions, and thus, in the case of heavy 
flavours, is applicable as long as the mass of the heavy flavour can be considered sufficiently large. The 
perturbative QCD cross section for heavy flavour production has been computed to next-to-leading order 
accuracy (i.e. 0{afj) a long time ago ^ [0] [§] [Q] [§] [§] [0] [§ [g] [[k| [juj], and a large amount of 
experimental and theoretical work has been done in this field. A relatively recent account of the status of 
this field can be found in ref. [12]. It can be said that qualitatively the QCD description of heavy flavour 
production seems to be adequate also for charm production, while quantitatively large uncertainties are 
present in the calculation of the charm and bottom cross section. Only for a quark as heavy as the top 
quark the perturbative calculation seems (up to now) to predict the cross section with a good accuracy. 
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Table 1 : Dependence of the b cross section on scale choices. 
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Large uncertainties are also found in the calculation of the bottom production cross section at the 
LHC. The largest uncertainty is due to unknown higher order effects, and it is traditionally quantified by 
estimating the scale dependence of the cross section when the renormalization and factorization scales 
are varied by a factor of 2 above and below their central value, which is usually taken equal to the 
heavy quark mass. Since this uncertainty is due to a limitation in our current theoretical knowledge, it is 
hard to overcome. Other sources of uncertainty are related to theoretical and experimental errors in the 
parameters that enter the perturbative calculation: the value of the strong coupling constant, the heavy 
quark mass, and the parton density functions. 

We present here a benchmark study of b total cross sections at the LHC, using the FMNR package 
for heavy flavour cross sections ^ [JsJ] (the code for this package is available upon request to the authors). 
In the study we consider 

• The dependence of the total cross section on the choice of the factorization and renormalization 
scales. We will use the values fi = nib, 2m&, m&/2. 

• The dependence on the parton density parametrization. We will use the sets MRST Jl3|], MRST(<?t), 
MRST(p|), MRST(a s |) and MRST(a s t). The first set is used as reference set. MRST(#t) and 
MRST(<7J.) have extreme gluon densities, MRST(a s |)-MRST(a s T) have extreme values of the 
strong coupling constant: A 5 = 220 MeV for MRST, 164 MeV for MRST(a s |), 280 MeV for 
MRST(a s T). Cross section values obtained with the CTEQ4 [14] set are very similar to the MRST 
set. We have preferred to use the MRST sets because they gave us the possibility to perform a 
study of sensitivity to A and to variations in the gluon density. 

• The dependence on the b quark mass: mj, = 4.75 ± 0.25 GeV. 

Factorization and renormalization scale dependence of the total cross section at V~S = 14 TeV is reported 
in table |J where we have used the MRST parton densities, with A5 = 220 MeV, and we have fixed the 
b mass at the value m& = 4.75 GeV. Notice that: 

• If we keep /ip = MR' the full cross section variation is small (467 to 512 fib). 

• The largest cross section corresponds to large /ip and small fi^ 

• The smallest cross section corresponds to small /ip and large /xr 

This is understood since, at small x, the gluon density g(x) grows with the scale, and a s decreases with 
the scale. 

The dependence on the choice of parton density parametrization is shown in table ||. As one can 
see, the sensitivity to the variation of the gluon density is small. Apparently, the constraints from HERA 
data are strong enough in the x region where most of the b production takes place. The dependence upon 
the strong coupling constant is instead larger, and can increase the upper limit of the cross section by 
about 40%. The last two sets have A5 = 164 and 288 MeV respectively, corresponding to a s (Mz) = 
0.1125 and 0.1225, which is a reasonably large range. 



Table 2: Parton density dependence of total cross sections (in [ib). 
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Mass uncertainties are quite important, especially if mj, is allowed to take very small values. This 
can be seen from table |3[ We see that lowering the b mass from 4.5 down to 4 GeV raises the upper limit 



Table 3: Mass dependence of total cross sections (in [ib). 
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central 


lowest 
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of the cross section by about 50%. It is however unlikely that such small values are viable. A rough 
view of the status of the bottom mass determination is given in fig. ^[ which we obtained by taking the 
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Fig. 2: Different determinations of the b quark pole mass. 



various determinations of the MS bottom mass from the Particle Data Book, and rescale them by a factor 
of (1+0.09+0.06), to account for the two-loop correction needed to translate the MS mass into the pole 
mass. As one can see, not all determinations are consistent among each other. A critical review of all 
determinations is beyond the scope of this workshop. We should however point out that recent progress 
has been made in the bottom mass determination. The reader can find a summary of these issues and 
further references in ref. [p^]. It is argued there that the bottom mass is determined with higher precision 
in processes where it is probed at short distances, like in the T mass, or via sum rules applied to the 
bottom vector current spectral function in e + e~ annihilation. The mass extracted in this way can be 



reliably related to the so called MS mass. The relation of the MS mass to the pole mass is instead not so 



precise, because the perturbative expansion that relates the two quantities is not convergent. In ref. [ ]15[ ] 
a preferred value of fht,(mb) = 4.23 ± 0.08 is given, where rhb(mb) is the MS bottom mass at the scale 
of the bottom mass itself. The corresponding pole mass, obtained using the newly computed 3-loop 
relation between the MS and pole mass [ |i~6[ ] JT7|], is 4.98 ± 0.09 GeV. If one wanted to account for the 
uncertainties due to the lack of convergence of the perturbative expansion, the range obtained in this way 
should be enlarged by some amount, of the order of 100 MeV. The question arises whether it would be 
possible to eliminate this uncertainty by expressing the hadroproduction cross section in terms of the MS 
mass. In our view, the answer is most likely no, since the bottom hadroproduction cross section does not 
have the same inclusive character as the sum rules applied to the e + e _ bottom spectral function. 

In the present work we thus used the traditional range 4.5 GeV < < 5 GeV for the bottom 
pole mass in the hadroproduction process, keeping in mind that recent determinations seem to favour the 
upper region of this range. The sensitivity of the cross section to the bottom mass in this range is at most 
of ±10%, and it becomes much smaller if transverse momentum cuts are applied. Thus, as far as the 
LHC is concerned, this uncertainty is not very important. 

The largest uncertainy in the cross section comes from the scale uncertainty, which is a (rather 
arbitrary) method to assess the possible impact of unknown higher order corrections. In the following 
we report a brief discussion of the origin of these large corrections. Radiative corrections for the total 
cross section are usually parametrized as follows. The total cross sections aij for the various parton 
subprocesses (qq, qg, gg) have a perturbative expansions given by 
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where p = 4ml /s an ^ * * s me squared partonic center-of-mass energy. The functions /£• for the 
and gg subprocesses are displayed in fig. ||. Notice the behaviour near threshold 
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Fig. 3: Partonic cross section for the qq and gg subprocesses. 
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due to Coulomb 1/(3 singularities and to Sudakov double logarithms. Near threshold, these terms may 
require special treatement, such as resummation to all orders. Notice also the constant asymptotic be- 
haviour of ffj , which may cause problems far above threshold. 

Plotting the cross section as a function of the partonic energy s may help to understand the origin 
of large corrections. We find that radiative corrections are large near the production threshold. This 



a 

• i— i 



60 



50 



40 



30 



1 



b 20 



10 



1 1 I I I 1 1 1| 1 1 I I I 1 1 1 1 1 1 I I I 1 1 1| 1 1 TTTTfrf 



I TTTTTTj 1 — I I I I Mil 



Lower: Born 
Upper: Full 




10° 10 1 



I I 



s/4m 2 

Fig. 4: Partonic cross section as a function of 3 /4m 2 



10 5 10 6 



problem becomes more and more severe as we approach the production threshold. Thus, it is more 
important for production of b at fixed target energies, or for production of ti pairs at colliders. Techniques 
to resum these large corrections to all orders of perturbation theory, at the NLO level are available [18], 
but it is found that little improvement is achieved for the bottom cross section at collider energies. Large 
corrections are also found far above threshold. This effect is bound to become more and more pronounced 
in the high energy limit. In order to reduce the scale uncertainties coming from these corrections, one 
should resum them at the next-to-leading order level. This problem has been discussed in the literature, 
so far, only at the leading logarithmic level [19] [ ]20| ] [ p| ] P2Q. At the time of the completion of this 
workshop, no further progress has been achieved in this field. 

In fig. |5| we present a study of the scale dependence of the total cross section as a function of 
p. We find a large scale dependence near threshold, due to both renormalization and factorization scale 
variation, and a large scale dependence far from threshold. Here, the renormalization scale dependence 
plays a dominant role. Renormalization scale variations are mainly due to the large variation of the 
coupling constant in the C(a|) terms. Where radiative corrections are small, a reasonable scale compen- 
sation takes place. Thus, both the threshold and the high energy regions, where corrections are large, are 
strongly affected. Factorization scale variation has a strong impact on threshold corrections, while in the 
high energy region we observe some compensation. In fact, the cross section near threshold increases 
with /ip near threshold, while above threshold the p-p = m value is above both the fip = m/2 and the 
/iF = 2m curves, indicating the presence of some sort of compensation. As of now, it appears there- 
fore that a better understanding of the high energy region will not strongly reduce the scale uncertainty, 
although it might, of course, improve our confidence in the error band we quote. 

The study given here deals with total cross sections. It should be repeated with appropriate rapidity 
cuts, since this may reduce large effects due to the high energy limit. In general, we may expect that the 




Fig. 5: Scale dependence of the total cross section. 
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Fig. 6: Cross sections as in table |4| The four bars in each group correspond to different choices of e and (fc-r), in the order 
reported in the figure. 



cross section with rapidity and transverse momentum cuts may have smaller error bars than the total. 
It is particularly interesting to investigate directly cross sections for muons originating from B decays, 
since muons are often used as trigger objects for B physics. We have performed this study using a 
simple implementation of the B semileptonic decay in the FMNR program, that will be described in 
more details in the following subsections. The results are shown in table ^. The same results are also 
reported in fig. ^, since several features become more apparent there. First of all, we point out that, as 
expected, there is a considerable reduction in the scale dependence in these muon rates. This is mostly 
due to the presence of cuts in the transverse momentum of the muon, that increases the total transverse 



Table 4: Cross sections (in fib) for b— »/tt + X production, with a muon, or both, satisfying appropriate cuts. Only muons 
coming directly from B decays are included here. The calculation was performed using the CTQ4M partem densities. The 
upper number are the maximum, and the lower number the minimum of the values obtained by varying the scales in the usual 
way. The corresponding total cross sectios are 165 to 864 fih The B— »/i branching fraction was taken equal to 10.5%. Different 
values for the e parameter of the Peterson fragmentation function are assumed. The last two column show the impact of a rather 
large intrinsic transverse momentum of the incoming partons. 
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A: bb^n(\r]\ < 2.4, p T > 6) 


3.3 
1.06 


2.41 
0.81 


2.12 
0.72 


3.4 
1.06 


2.91 
0.94 


B: bb—>n(pr > 6) h(j>t > 3) 


0.76 
0.304 


0.52 
0.219 


0.45 
0.19 


0.67 
0.252 


0.54 
0.214 


C: bb^fi(p T > 6) e(p T > 2) 


1.18 
0.43 


0.83 
0.32 


0.71 
0.277 


1.1 

0.38 


0.92 
0.33 


D: bb^fi( PT > 7, \rj\ < 2.4) 


2.26 
0.78 


1.62 
0.58 


1.41 
0.5 


2.23 
0.73 


1.9 
0.63 


E: bb^niPT > 7, \rj\ < 2.4) 

fi(p T > 4.5, < \rj\ < 1.5) 


0.0304 
0.0146 


0.0203 
0.0102 


0.0174 
0.0087 


0.0232 
0.0105 


0.0188 
0.009 


F: bb^n(p T > 7, \rj\ < 2.4) 

H{PT > 3.6, 1.5 < M < 2) 


0.0101 
0.0045 


0.0075 
0.0032 


0.0068 
0.0026 


0.0096 
0.0035 


0.0076 
0.00281 


G: bb^ii{p T > 7, \q\ < 2.4) 

fi{p T > 2.6, 2 < \r)\ < 2.4) 


0.0105 
0.0038 


0.0073 
0.00263 


0.0053 
0.00219 


0.0082 
0.00251 


0.0062 
0.0024 


H: bb^fi( PT > 1, 2 < \r)\ < 6) 


19.3 
5.4 


18.8 
5.3 


18.6 
5.2 


19.1 
5.4 


18.9 
5.3 


I: bb^fi(p T > 2, 2 < |r/| < 6) 


10.2 
2.94 


9.1 
2.65 


8.6 
2.51 


10.6 
3.11 


10. 
2.96 



energy that characterizes the cross section. Thus, while the ratio of the upper to the lower limit of the 
cross section is above a factor of 5 in the total rate, it is between a factor of 2 and 4 in the muon rates. 
The smallest values are achieved for the highest momentum cuts. A non perturbative fragmentation 
function of the Peterson form was also included in the calculation, with e parameter taking the values 
(i.e., no fragmentation), 0.002 and 0.006. More details on its implementation are given in the following 
subsections. Observe that for softer fragmentation functions (i.e. larger e parameter) the uncertainty is 
reduced, since they imply higher quark momenta. The reduction in the scale uncertainty is obtained at 
the price of introducing a sensitivity to the fragmentation function parameter. We considered as realistic 
values of e between 0.002 and 0.006. The corresponding variation of the cross section is not large. The 
impact of an intrinsic transverse momentum of the incoming partons (see the following subsections) is 
also studied. We have chosen the unrealistically large value (kx) = 4 GeV just to show that its effect is 
in all cases not a dramatic one. 

2.2 Transverse momentum spectrum 

2.21 Benchmark single-inclusive distributions 

The fixed-order, NLO result for single-inclusive b production has several limitations in different regions 
of the phase space. In particular, one should be aware of the high-energy limit problem when px is 
small compared to the incoming energy, of the logarithms of nib/pT for high transverse momenta, and 
of further problems when approaching the threshold region. All these issues will be discussed in some 
detail in the next Sections. However, the fixed-order calculation at NLO provides a useful starting point 
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Fig. 7: Differential cross section for heavy flavour production vs. pr, for different rapidities. 
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Fig. 8: Differential cross section for heavy flavour production vs. y for different pr values, as given in tables 



for estimating the differential cross section. At this time, it is probably not useful to perform a cross 
section study with different sets of parton densities, and for different values of the b mass. We limit 
ourselves to the MRST set, and we only study the scale dependence of the cross section. We do not 
include, at this stage, fragmentation effects, which, as shown in the following Sections, can be easily 
accounted for. In tables ||-|8| we collect the results of this study. The central values we obtained are also 
plotted in figs. ^ and [| so that the wide kinematic range of heavy flavour production can be appreciated 
by a glance. More detailed rapidity distributions at low momenta are shown in fig. |9[ First of all, we see 
that the differential cross section spans many orders of magnitude. At a luminosity of 10 34 cm _2 sec _1 
each fib of cross section corresponds to 10 4 events per second, or (roughly) 10 11 events per year. Thus, 
at the level of 10 -11 in the plot there should be one event per year per bin of px and y. The px spectrum 
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Fig. 9: Differential cross section for heavy flavour production vs. y for different pr values. 



starts to drop fast for pt larger than the heavy quark mass, dropping even faster as the threshold region 
is approached. The rapidity distributions have the typical shape of a wide plateau, dropping at the edge 
of the phase space, and becoming narrower for larger transverse momenta. At the LHC the gluon fusion 



production mechanism is dominant, as can bee seen in fig. 10. There one can see that the quark- antiquark 
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Fig. 10: The total, qq and qg components of the differential cross section for heavy flavour production. 



annihilation component is below the gluon fusion component by more than one order of magnitude in the 
Pt range considered, while the quark-gluon term becomes more important at larger pr- We remind the 
reader that the cross section for qq—>b + X is not included in the NLO calculation. One may thus worry 
about a loss of accuracy in the result, since the quark-quark luminosity at the LHC are by far the largest 
for high transverse momenta b production. This problem, however, is dealt with appropriately in the 



resummation formalism for high pr heavy flavour production, where a quark-quark fusion contribution 
does indeed appear. 



Table 5: Differential cross section da/dp^dy for single inclusive b production at the LHC, forpr from 1 to 80 GeV and y from 
to 4. The table was computed with the MRST parton density set, for m;, = 4.75 GeV. The central value was obtained with 
the factorization and renormalization scale set to \Jp\-\~ m 2 . The upper and lower values give the maximal variation when 
varying the scales independently by a factor of 2 above and below the central value. Cross sections are in fib; each element in 
a row should be multiplied by the common scale factor in the left column. 
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34.1 


31 


26.1 


PT = 


2 


6.7 


6.6 


6.5 


6.3 


6 


5.7 


5.2 


4.56 


10" 


x x 


19.2 


18.8 


18.2 


17.4 


16.4 


15.1 


13.5 


11.5 






30.3 


29.7 


29 


28 


26.7 


24.9 


22.6 


19.2 


PT = 


3 


5.4 


5.3 


5.2 


5 


4.76 


4.43 


4.01 


3.48 


10" 


1 x 


14 


13.6 


13.2 


12.6 


11.8 


10.8 


9.5 


8 






21.8 


21.3 


20.7 


19.9 


18.8 


17.4 


15.6 


13.2 


PT = 


4 


4.04 


3.96 


3.85 


3.7 


3.48 


3.21 


2.87 


2.45 


io- 


x x 


9.9 


9.6 


9.3 


8.8 


8.2 


7.4 


6.5 


5.4 






15.9 


15.5 


15 


14.2 


13.3 


12.1 


10.6 


8.9 


PT = 


5 


2.91 


2.84 


2.75 


2.62 


2.45 


2.24 


1.97 


1.64 


10" 


1 x 


6.6 


6.4 


6.2 


5.8 


5.3 


4.79 


4.13 


3.37 






11.1 


10.7 


10.3 


9.7 


9 


8 


6.9 


5.7 


PT = 


7 


14 


13.6 


13.1 


12.3 


11.3 


10.1 


8.6 


6.9 


10" 


2 x 


29.4 


28.2 


26.9 


25 


22.7 


19.8 


16.6 


13 






51 


48.5 


46 


42.7 


38.5 


33.5 


27.9 


21.7 


PT = 


10 


4.72 


4.54 


4.31 


3.99 


3.58 


3.09 


2.52 


1.9 


10" 


2 x 


9.1 


8.7 


8.2 


7.5 


6.6 


5.6 


4.49 


3.32 






15.6 


14.8 


13.9 


12.6 


11.1 


9.3 


7.4 


5.4 


PT = 


15 


10.3 


9.8 


9.1 


8.2 


7.1 


5.9 


4.49 


3.11 


10" 


3 x 


17.8 


16.8 


15.6 


14 


12 


9.7 


7.3 


4.91 






29.3 


27.5 


25.4 


22.6 


19.2 


15.4 


11.4 


7.6 


PT = 


20 


2.99 


2.81 


2.6 


2.3 


1.94 


1.53 


1.1 


0.7 


10" 


3 x 


4.86 


4.55 


4.16 


3.65 


3.04 


2.36 


1.67 


1.04 






7.7 


7.1 


6.5 


5.7 


4.69 


3.61 


2.52 


1.55 


PT = 


30 


4.53 


4.2 


3.79 


3.25 


2.61 


1.91 


1.25 


0.68 


10" 


4 x 


6.8 


6.3 


5.6 


4.77 


3.78 


2.74 


1.75 


0.93 






10.1 


9.3 


8.3 


7 


5.5 


3.94 


2.49 


1.3 


PT = 


40 


10.9 


10 


8.8 


7.4 


5.7 


3.9 


2.31 


1.08 


10" 


5 x 


15.5 


14.1 


12.5 


10.3 


7.8 


5.3 


3.1 


1.41 






22.3 


20.2 


17.8 


14.6 


11 


7.4 


4.23 


1.92 


PT = 


50 


3.47 


3.13 


2.72 


2.21 


1.63 


1.06 


0.57 


0.218 


10" 


5 x 


4.77 


4.28 


3.71 


2.99 


2.18 


1.4 


0.74 


0.287 






6.7 


6 


5.2 


4.12 


2.98 


1.89 


1.01 


0.395 


pr = 


60 


13.2 


11.8 


10.1 


8 


5.7 


3.48 


1.66 


0.53 


10" 


6 x 


17.8 


15.8 


13.5 


10.6 


7.5 


4.52 


2.18 


0.71 






24.3 


21.5 


18.3 


14.3 


10 


6.1 


2.97 


0.98 


PT = 


80 


2.77 


2.42 


2.01 


1.5 


0.98 


0.53 


0.207 


0.0449 


io- 


6 x 


3.59 


3.12 


2.59 


1.95 


1.28 


0.69 


0.273 


0.06 






4.76 


4.14 


3.45 


2.6 


1.72 


0.94 


0.376 


0.083 



Table 6: As in table pi for y from 4.5 to 6.5. 



y 


A C 

4.5 


5 


5.5 


6 


6.5 


PT = 1 


4.51 


3.46 


2.36 


1.36 


0.59 


10 x x 


11.6 


8.4 


5.3 


2.87 


1.15 




20 


14.2 


8.5 


4.47 


1.8 


p T = 2 


3.74 


2.82 


1.88 


1.04 


0.431 


10 L x 


9 


6.4 


4.04 


2.1 


0.8 




1 A C 

14.b 


10.2 


0.3 


3.2 


1.23 


PT = 3 


2.8 


2.06 


1.32 


0.69 


0.261 


10 A x 


6.2 


4.35 


2.65 


1.31 


0.458 




10.1 


fa.8 


4.07 


1.97 


0.7 


PT — 4 


1.93 


1.37 


0.84 


0.415 


0.137 


10 x x 


4.08 


2.78 


1.62 


0.75 


0.231 




b.7 


A A A 

4.44 


2.53 


1.15 


a or 

0.35 


Pr = 5 


1.27 


0.87 


0.51 


0.233 


0.065 


10 x x 


2.52 


1.67 


0.93 


0.404 


0.105 




4.19 


2.71 


1.5 


O.o3 


O.lo 


p T = 7 


5.1 


3.24 


1.71 


0.65 


0.12 


10 2 x 


9.3 


5.8 


2.93 


1.06 


0.178 




1 r J 

15.4 


9.4 


A PC 

4.00 


l.o3 


0.268 


= 10 


1.29 


0.75 


0.333 


0.09 


6.4 10" 3 


10 


2.2 


1.23 


0.53 


0.134 


8.8 10~ 3 




o r a 

3.54 


1.95 


0.81 


0.2 


0.0128 


= 15 


1.86 


0.88 


0.279 


0.0324 


1.2 10 -4 


10 3 x 


2.86 


1.31 


0.395 


0.0431 


1.57 10~ 4 




A OP 

4. 3o 


1.95 


0.57 


O.Ool 


2.25 10 


= 20 


0.371 


0.143 


0.0288 


8.8 10" 4 


1.17 10~ 12 


10 3 x 


0.54 


0.199 


0.0379 


1.18 10~ 3 


a A -t r\—i 9 

1.4410 12 




0.78 


n 001 

0.281 


0.053 


1 cc\ 1 n — 3 

l.o910 


1.57 10 


PT = 30 


0.278 


0.062 


3.28 10~ 3 


1.14 10~ 7 





r. .A 

10 4 x 


0.369 


0.082 


4.39 10~ 3 


1.55 10~ 7 







0.5 


0.114 


0.2 10 


2.02 10 





PT = 40 


0.316 


0.036 


2.48 10~ 4 








10 5 x 


0.418 


0.0481 


3.24 10 -4 










0.58 


O.Otw 


A A Q i n — 4 
4.43 10 








Pt = 50 


0.0457 


2.25 10" 3 


6.4 10" 8 








10" 5 x 


0.061 


3.02 10" 3 


8.2 10~ 8 










0.085 


4.2 10" 3 


1.07 10" 7 








Pt = 60 


0.076 


1.23 10" 3 











10~ 6 x 


0.102 


1.64 10~ 3 













0.143 


2.25 10" 3 











Pt = 80 


2.48 10" 3 


2.93 10" 7 











10" 6 x 


3.29 10~ 3 


3.9 10~ 7 













4.56 10" 3 


4.78 10" 7 












Table 7: As in table |, for p T from 100 to 300. 



y 





1 


1.5 


2 


2.5 


3 


3.5 


4 


4.5 


p T = 100 


7.8 


6.6 


5.3 


3.82 


2.33 


1.11 


0.354 


0.0484 


6.9 10" 4 


10~ 7 x 


10 


8.5 


6.9 


4.97 


3.05 


1.47 


0.469 


0.064 


910" 4 




13.3 


11.4 


9.3 


6.7 


4.14 


2.02 


0.65 


0.09 


1.23 10" 3 


p T = 120 


2.66 


2.22 


1.75 


1.2 


0.69 


0.291 


0.073 


5.8 10" 3 


9.510"^ 


10~ 7 x 


3.44 


2.89 


2.28 


1.57 


0.9 


0.387 


0.098 


7.7 10" 3 


1.11 10" 5 




4.59 


3.87 


3.06 


2.13 


1.23 


0.53 


0.135 


0.0106 


1.54 10" 5 


p T = 140 


10.5 


8.7 


6.7 


4.41 


2.35 


0.88 


0.172 


7 10 -3 


4.61 10" 8 


10~ 8 x 


13.7 


11.3 


8.7 


5.8 


3.1 


1.17 


0.229 


9.1 10~ 3 


5.4 10~ 8 




18.4 


15.2 


11.8 


7.9 


4.24 


1.62 


0.317 


0.0125 


6.1 10" 8 


p T = 160 


4.68 


3.79 


2.84 


1.81 


0.9 


0.296 


0.0433 


7.6 10" 4 





10" 8 x 


6.1 


4.96 


3.73 


2.38 


1.19 


0.396 


0.058 


9.9 10" 4 







8.2 


6.7 


5.1 


3.24 


1.63 


0.55 


0.08 


1.33 10- 3 





p T = 180 


22.6 


18 


13.2 


8.1 


3.74 


1.08 


0.115 


6.8 10~ 4 





10~ 9 x 


29.5 


23.6 


17.4 


10.7 


4.96 


1.43 


0.152 


8.7 10~ 4 







39.8 


31.9 


23.6 


14.6 


6.8 


1.99 


0.21 


1.13 10 -3 





p T = 200 


11.7 


9.2 


6.6 


3.87 


1.67 


0.415 


0.0312 


3.29 10" 5 





10~ 9 x 


15.3 


12.1 


8.7 


5.1 


2.21 


0.55 


0.0413 


4.68 10" 5 







20.6 


16.3 


11.8 


7 


3.03 


0.76 


0.056 


5.4 10" 5 





p T = 220 


6.4 


4.93 


3.46 


1.95 


0.78 


0.167 


8.4 10~ 3 


4.96 10- y 





10~ 9 x 


8.4 


6.5 


4.58 


2.59 


1.04 


0.223 


0.0111 


6.5 10~ 7 







11.3 


8.8 


6.2 


3.54 


1.43 


0.307 


0.015 


7.1 10" 7 





p T = 240 


3.66 


2.79 


1.91 


1.04 


0.385 


0.07 


2.24 10" 3 


2.8210~ 1U 





10~ 9 x 


4.8 


3.67 


2.52 


1.37 


0.51 


0.092 


2.92 10" 3 


810- 10 







6.5 


4.99 


3.45 


1.88 


0.7 


0.127 


3.91 10" 3 


9.7 10~ 10 





p T = 260 


21.8 


16.3 


11 


5.7 


1.96 


0.298 


5.7 10" 3 








10- 10 x 


28.7 


21.6 


14.5 


7.6 


2.61 


0.396 


7.4 10~ 3 










38.8 


29.3 


19.8 


10.3 


3.58 


0.54 


9.6 10" 3 








p T = 280 


13.4 


9.9 


6.5 


3.24 


1.02 


0.13 


1.32 10" 3 








10- 10 x 


17.7 


13.1 


8.6 


4.32 


1.37 


0.172 


1.73 10- 3 










23.9 


17.8 


11.8 


5.9 


1.87 


0.234 


2.18 10" 3 








p T = 300 


8.5 


6.2 


3.98 


1.9 


0.55 


0.057 


2.7 10" 4 








10- 10 x 


11.2 


8.2 


5.3 


2.52 


0.73 


0.076 


3.67 10" 4 










15.2 


11.1 


7.2 


3.46 


1.01 


0.103 


4.4 10" 4 









Table 8: As in tabled forpT from 320 to 500. In the entries with a * the cross section is too small to be computed reliably. 



y 





1 


1.5 


2 


2.5 


3 


3.5 


PT = 


320 


5.5 


3.97 


2.5 


1.14 


0.304 


0.0256 


4.36 10 -5 


10" 


i n 

iu x 


7.3 


5.2 


3.3 


1.51 


0.403 


0.0334 


6.5 10" 5 






9.9 


7.2 


4.52 


2.08 


0.55 


0.0453 


7.3 10~ 5 


PT = 


340 


3.67 


2.6 


1.6 


0.7 


0.17 


0.0114 


5.1 10 _e 


icr 


i n 

iu x 


4.85 


3.44 


2.12 


0.93 


0.226 


0.0149 


7.2 10" 6 






6.6 


4.69 


2.89 


1.27 


0.309 


0.0199 


7.710 


PT = 


360 


2.49 


1.74 


1.04 


0.437 


0.097 


5.1 10" 3 


3.56 10~ 7 


icr 


i n 

iu x 


3.29 


2.31 


1.38 


0.58 


0.128 


6.6 10~ 3 


4.33 10 -7 






4A7 


3.13 


1.89 


0.79 


0.175 


8.7 10 -3 


4.63 10 -7 


PT = 


380 


17.1 


11.8 


6.9 


2.78 


0.56 


0.0223 


6.410~ s 


10" 


11 X 


22.8 


15.7 


9.2 


3.7 


0.74 


0.0288 


1.24 10 -7 






30.9 


21.4 


12.6 


5.1 


1 


0.038 


1.46 10 


PT = 


400 


12 


8.2 


4.69 


1.79 


0.326 


9.7 10" 3 




10" 


11 X 


16 


10.8 


6.2 


2.38 


0.43 


0.0124 


* 






2L.7 


14.8 


8.5 


3.2b 


0.58 


U.U163 




PT = 


420 


8.6 


5.8 


3.21 


1.17 


0.192 


4.16 10" 3 




lO - 


11 X 


11.3 


7.6 


4.25 


1.55 


0.251 


5.2 10~ 3 


* 






15.5 


10.4 


5.8 


2.13 


0.342 


6.8 10" 3 




PT = 


440 


6.2 


4.09 


2.23 


0.78 


0.114 


1.73 10" 3 





10" 


11 X 


8.2 


5.4 


2.96 


1.03 


0.149 


2.14 10" 3 









11.2 


7.4 


4.04 


1.41 


0.201 


2.79 10" 3 





PT = 


460 


4.53 


2.94 


1.57 


0.52 


0.068 


6.9 10~ 4 





10" 


11 X 


6 


3.9 


2.08 


0.69 


0.089 


8.5 10~ 4 









8.2 


5.3 


2.84 


0.94 


0.119 


1.09 10" 3 





PT = 


480 


3.34 


2.14 


1.11 


0.351 


0.0407 


2.61 10" 4 





10" 


11 X 


4.43 


2.84 


1.48 


0.465 


0.053 


3.16 10" 4 









6 


3.87 


2.01 


0.63 


0.071 


4.07 10" 4 





PT = 


500 


2.49 


1.57 


0.8 


0.239 


0.0244 


9.3 10" 5 





10" 


11 X 


3.31 


2.09 


1.06 


0.316 


0.0319 


1.1 HT 4 









4.51 


2.85 


1.44 


0.429 


0.0423 


1.45 10" 4 






2.22 Understanding Tevatron data 

It is well known that Tevatron data for the integrated transverse momentum spectrum in b production are 
systematically larger than QCD predictions. This problem has been around for a long time, although it 



has become less severe with time. The present status of this issue is summarized in fig. 11. A similar 
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Fig. 1 1 : The integrated pr distribution for single b production measured at the Tevatron, and the corresponding QCD prediction. 



discrepancy is also observed in UA1 data (see ref. [12] for details). 

The theoretical prediction has a considerable uncertainty, which is mainly due to neglected higher- 
order terms in the perturbative expansion. In our opinion, it is not unlikely that we may have to live 
with this discrepancy, which is certainly disturbing, but not strong enough to question the validity of 
perturbative QCD calculations. In other words, the QCD 0(a\) corrections for this process are above 
100% of the Born term, and thus it is not impossible that higher order terms may give contributions of 
the same size. Nevertheless, it is useful to look for higher-order perturbative effects and non-perturbative 
effects that may enhance the cross section. 

For values of pr much larger than the b quark mass, large logarithms of the ratio pr / mj, arise in the 
coefficients of the perturbative expansion. Techniques are available to resum this class of logarithms to all 



orders. In ref. [ |23| ] the NLO cross section for the production of a massless parton i (a gluon or a massless 
quark) has been folded with the NLO fragmentation function for the transition i—>b [24]. The evolution 



of the fragmentation functions resums all terms of order a" log™ (pr / 'mi,) and a n log n (pr/mi,). All 
the dependence on the 6-quark mass lies in the boundary conditions for the fragmentation function. The 
result is then matched with the full NLO cross section, which contains the exact dependence on m^ up 
to order af , in a way that avoids double counting. Corrections to the result of ref. j f23| ] are either of order 

a\ log* (pr I mf,) , with i < 2, or of order a| times positive powers of m&/ \/Pr + m 1- 



Figures |12j-|13| show the differential and integrated 6-quark pr distribution obtained in the frag- 
mentation function approach of ref. [E3p, compared to the standard fixed-order NLO result. It should be 
noted that for high transverse momenta the scale dependence is significantly reduced with respect to the 



fixed-order calculation. Furthermore, it can be seen from fig. |f| that, for 10 J> px iJ 30 GeV, the result 
of the fragmentation-function approach lies slightly above the fixed-order NLO calculation. This has 



been interpreted in ref. [23] as an evidence for large, positive higher order corrections. Unfortunately, 
their effect is not easy to quantify. These higher order terms are in fact computed in a massless approx- 



imation, and thus fail at low transverse momenta. In figs. 12 13 these terms are suppressed by a factor 
that becomes smaller and smaller at low p T . A more detailed discussion of this point can be found in 



10 c 



10- 



% 10 4 

o 

& 10 J 

£ io 1 



T3 10 







10" 



10 



-2 



— i — i — i — I — r~ 




QCN NLO 

QCD NLO + resummation 



E CM =1.8 TeV, CTEQ3M 
m b = 5 GeV, y=0 



40 60 
Pt (GeV) 



100 



Fig. 12: Single-inclusive pr distribution for b production at the Tevatron energy: pure QCD and resummed results. 
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Fig. 13: Integrated pr distribution for b production at the Tevatron energy: pure QCD and resummed results. 



the original reference. Here, we simply conclude that some evidence for large higher order terms in the 
intermediate transverse momentum region is present, although difficult to quantify. 

Finally, notice that the overall effect of the inclusion of higher-order logarithms is a steepening 
of the pt spectrum. This is quite natural, since multiple radiation is accounted for in the resummation 
procedure. 



It has been argued that an intrinsic transverse momentum for the incoming partons may explain the 
discrepancy observed at the Tevatron. In fact, large values (up to 4 GeV) of the average transverse mo- 
mentum of the incoming partons have been invoked to explain direct photon production data J25|]. Such 
large values, much larger than typical QCD scales, are clearly incompatible with the usual application of 
perturbative QCD. Thus, evidence for such a large intrinsic transverse momentum cannot be claimed on 
the basis of a single observable. In other words, we would need evidence from several observables, all 
leading to a similar value of the intrinsic kr, before we accept such a flaw in the usual perturbative QCD 
description. Nevertheless, in the following we will perform the exercise of applying very large intrinsic 
transverse momenta to the heavy flavour production process. This procedure will lead to an increase in 
the b transverse momentum spectrum. We will also show, however, that other variables, very sensitive to 
an intrinsic transverse momentum, that should be strongly affected, do not show any evidence of that. 

There are several possible ways to implement the presence of a non-zero transverse momentum of 
the colliding partons, and the choice is, to a large extent, arbitrary. We implemented it in the FMNR code 
in the following way. We call p T (QQ) the total transverse momentum of the pair. For each event, in 
the longitudinal centre-of-mass frame of the heavy-quark pair, we boost the QQ system to rest. We then 
perform a transverse boost, which gives the pair a transverse momentum equal to p T (QQ) + fc T (l) + 
k T (2); k T (l) and /c T (2) are the transverse momenta of the incoming partons, which are chosen randomly, 
with their moduli distributed according to 



1 



1 dN 



N dk% {kl) 



exp(- 



-k\j (k^)). 



(2) 



The reader can find more details in ref. [ 12] 



In fig. |14| we show the effect of an intrinsic kx generated in this way, with the (unphysically large) 
choice (k T ) = 4 GeV (in fig. 14, the sensitivity to the parameter in the fragmentation function is also 
shown; fragmentation will be discused in more detail in the next subsection.) We see that, for p™ n < 20 
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Fig. 14: The b cross section at the Tevatron: the effect of a large intrinsic transverse momentum, and the sensitivity to the 
fragmentation parameter e^. 



GeV, the k? effect is sizeable, even in the presence of fragmentation, provided we allow for unphysically 
large intrinsic kr- 

It is fair to ask whether such large values are compatible with other observables. There is a particu- 
lar class of observables that are particularly sensitive to the intrinsic transverse momentum. One example 



is the azimuthal distance A</> between the directions of the produced b and b. The A(f> distribution is triv- 
ial at leading order: b and b are emitted back-to-back, and therefore 



da 

~dK4> 



OC 5{(j) — 71"). 



(3) 



An intrinsic kx of the colliding partons has the effect of smearing the S function. For (kx) = 4 GeV the 



effect is quite dramatic, as can be seen in fig. 15. Is such an important effect consistent with the observed 
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Fig. 15: The b-b azimuthal correlation at the Tevatron: the effect of a large intrinsic transverse momentum, and the sensitivity 
to the fragmentation parameter e b . 



azimuth correlations? The CDF and DO collaborations have measured the azimuthal correlation of muon 
pairs produced in b decays. In order to compare with these data sets, we have implemented in the FMNR 
code the semileptonic decay of b quarks. We have assumed that the muon energy is distributed according 



to the prediction of the spectator model [26] with massless leptons. We have also checked that the muon 



energy distribution given by Pythia leads to similar results. Our results are shown in figs. 16 and 17 
where CDF and DO data are superimposed to the perturbative QCD prediction, with and without an 
intrinsic kx with (kx) = 4 GeV. Tevatron data do not seem to favour such a large intrinsic transverse 
momentum. The measured distributions are more peaked at A<p = it than the theoretical curve with 
(kx) = 4 GeV. The effect of Peterson fragmentation is also shown in both cases. We thus conclude that 
the data does not seem to favour large kx effects. 

2.23 Single-incusive distributions and correlations at the LHC 

In this subsection, we will follow the pragmatic assumption that the discrepancy observed at the Tevatron 
may either be attributed to a problem in the overall normalization of the cross section, or to the presence 
of effects, either perturbative or not, that distort the spectrum. We will continue to model these effects 
as fragmentation effects and intrinsic transverse momentum effects, and see if the LHC can distinguish 



among the two. In fig. |18| we plot the b cross section with a transverse momentum cut. From the figure it 
is quite clear that the effects of fragmentation and the effects of an intrinsic transverse momentum kick 
manifest themselves in quite a different way. In particular, at p T > 20 GeV even the effect of a very 
large transverse momentum kick is small, while fragmentation has a strong impact. On the other hand, 
the transverse momentum kick increases the cross section in the intermediate p T region, with a maxi- 
mum around 7 GeV. The p T coverage offered by the combined LHC experiments will allow an effective 



10 



10 



10" 



-2 



-3 _ 



10 



-5 _ 




solid: NLO QCD 

lower dashed: <k T >=0, e b =0.002 
lower dotted: <k T >=0, E b =0.006 
upper dashed: <k T >=4 GeV, e b =0.002 
upper dotted: <k T >=4 GeV, € b =0.006 - 



Fig. 16: CDF results on azimuthal correlations compared with the perturbative calculation, with and without intrinsic kr. 
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Fig. 17: DO results on azimuthal correlations compared with the perturbative calculation, with and without intrinsic kr 



discrimination of the two kinds of effects. For completeness, we also show in fig. 19 a comparison 



of the fixed-order calculation of the single-inclusive spectrum, using the fixed-order calculation in two 
different schemes for the light flavour, and the matched-resummed result. As in the Tevatron case, the 
band obtained with the resummation procedure is much narrower at large transverse momentum. The 
corresponding uncertainty does not include fragmentation function uncertainties, that will be discussed 
in more detail further on. 

As an example of what could be discriminated at the LHC using correlations, we present in fig. |(] 
the azimuthal correlation of the muons coming from semileptonic B decays, using typical cuts that are 
implemented in the LHC experiments triggers for B studies. The curves are obtained with different 
values of the e& parameter for the fragmentation function, and with or without a very large intrinsic 




Fig. 18: Cross section with a transverse momentum cut at the LHC. 
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Fig. 19: Differential cross section for b production at the LHC. The bands are obtained by varying independently the renormal- 
ization and factorization scales by a factor of 2 above and below the central value, which is the b transverse mass. 



transverse momentum for the incoming partons. As one can expect, the 65 parameter affects only the 
total rate in this case, while the primordial transverse momentum has a considerable effect on the shape 
of the distribution. This example shows that, even with very simple experimental setup, at the LHC it 
will be possible to test important features of the differential distributions. 



2.3 Fragmentation function formalism 

In analogy with the case of charm production, the agreement between theory and data improves if one 
does not include any fragmentation effects. It is then natural to ask whether the fragmentation functions 
commonly used in these calculations are appropriate. Following the LEP measurements, fragmentation 
functions have appeared to be harder then previously thought. It will be interesting to see whether SLD 
new data [27] will help in clarifying this issue. 

The effect of a non-perturbative fragmentation function on the p-r spectrum is easily quantified if 
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Fig. 20: Azimuthal correlations in muon pairs at the LHC. 



one assumes a steeply-falling transverse momentum distribution for the produced b quark 

' la -W , (4) 



dpT 

The corresponding distribution for the hadron is 



dpi 



A J p? M 5(pt - zpt)D(z) dz dp T = Apr M dzz M ~ 1 D(z) . (5) 



We can see that the hadron spectrum is proportional to the quark spectrum times the M th moment of the 
fragmentation function D{z). Thus, the larger the moment, the larger the enhancement of the spectrum. 

In practice, the value of M will be slightly dependent upon pr- We thus define a px dependent M 

value 

d\oga{p T >p c ^) _ cut 



dlogp^ 1 

and 



-M{pf l ) + 1 (6) 



a had ( PT > PT) = °iPT > P c r nt ) x f dzz M ^)~l D ( z ) . ( 7 ) 

Jo 

This gives an excellent approximation to the effect of the fragmentation function, as can be seen from 



fig. 21 



Since the second moment of the fragmentation function is well constrained by e + e~ data, it is 
sensible to ask for what shapes of the fragmentation function, for fixed (z), one gets the highest value 
for (z M ~ 1 ). We convinced ourselves that the maximum is achieved by the functional form 

D{z) = A5{z) + B5{1 - z) (8) 

which gives 
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Fig. 21: The effect of a Peterson fragmentation function of the inclusive b cross section. The (red) dotted lines correspond to 
the approximation eq. ([]), and are almost indistinguishable from the exact results. 



This is however not very realistic: somehow, we expect a fragmentation function which is concentrated 
at high values of z, and has a tail at small z. We convinced ourselves that, if we impose the further 
constraint that D(z) should be monotonically increasing, one gets instead the functional form 

D{z) = A + B8{1 - z) , (10) 

which gives 

( z) = A / 2 + B . {z M-u = A/M + B 

{) A + B ' ' A + B ' U } 

We computed numerically the M th moments of the Peterson form, 

D(z) oc — ^ (12) 

z 1- * 



1-2 

of the form 

D(z) oc z a (l - zf (13) 

for f3 = 1 (Kartvelishvili), for which 

M-n _ r(a + M)T(a + (3 + 2) 
{ ' T{a + l)T(a + l3 + M + l) ' ^ ; 

of the form of Collins and Spiller 

D(z) oc ^ ; J —— 2 (15) 

and of the form in eq. (|I0|), at fixed values of (z) con - esponding to the choices e& = 0.002 and 0.006 
in the Peterson form. We found that the pr distribution at the Tevatron, for pt in the range 10 to 100 
GeV, behaves like p^ M , with M around 5. Therefore, we present in tables |9] and 10 values of the 4 th , 



Table 9: Values of the 4 th , 5 th and 6 th moment, at fixed {z} (corresponding to et = 0.002 in the Peterson form), for different 
forms of the fragmentation function. 



(z) = 0.879 


M = 4 


M = 5 


M = 6 


Peterson 


0.711 


0.649 


0.595 


Kartvelishvili 


0.694 


0.622 


0.562 


Collins-Spiller 


0.729 


0.677 


0.633 


Maximal (eq. (10))) 


0.818 


0.806 


0.798 



Table 10: Values of the 4 th , 5 th and 6 th moment, at fixed (z) (corresponding to ei, — 0.006 in the peterson form), for different 
forms of the fragmentation function. 



(z) = 0.828 


M = 4 


M = 5 


M = 6 


Peterson 


0.611 


0.535 


0.474 


Kartvelishvili 


0.594 


0.513 


0.447 


Collins-Spiller 


0.626 


0.559 


0.505 


Maximal (eq. (10))) 


0.742 


0.724 


0.713 



5 and 6 moments of the above-mentioned fragmentation functions. We thus find that keeping the 
second moment fixed the variation of the hadronic px distribution obtained by varying the shape of the 
fragmentation function among commonly used models is between 5% and 13% for both values of e^. It 
thus seems difficult to enhance the transverse momentum distribution by suitable choices of the form of 
the fragmentation function. With the extreme choice of eq. (|10|), one gets at most a variation of 50% for 
the largest value of ej, and M. It would be interesting to see if such an extreme choice is compatible with 
e + e~ fragmentation function measurements. 

3. A STUDY OF HEAVY QUARK NON-PERTURBATIVE FRAGMENTATION IN HERWIG| 

In this Section we present the results of a phenomenological study of the non-perturbative hadronization 
of 6-quarks. According to the standard QCD picture, distributions for an observable hadron H can be 
computed by convoluting the short-distance cross section a(p) with a fragmentation function D^\z) 
that describes the way in which the heavy quark h hadronizes into H: 

da H {p) = J dzD^\z)da{p/z). (16) 

The precise definition of Dg' (z) depends on how much of the heavy quark evolution after its production 
is absorbed into the perturbative part a(p), and how much is assigned to the non-perturbative component 
parameterised by D^\z). Since perturbation theory (PT) is well defined for a massive quark, the stan- 
dard prescription is to absorb into a(p) not only the hard matrix elements, but also the perturbative part 
of the fragmentation function, defined by the evolution in Q 2 down to a scale equal to the heavy quark 
mass m/j. Dw(z) will therefore account for the transition of an "on-shell" quark h into the hadron H. 
The assumptions built into eq. ( [I6| ) are that (z) depends neither on the type of hard process, nor on 

the scale at which h was produced. Under these assumptions, D^\z) can be extracted from data in one 
given reaction (typically, e + e~), and eventually used to predict the cross section in some other reaction 
(pp, DIS and so on). 

2 Section coordinators: S. Frixione and M.L. Mangano 



QCD factorization theorems indicate that this universality of (z) holds in the asymptotic limit, 
and up to corrections of order rn^/Q, Q being the scale of the hard process. The size of these corrections 
cannot be calculated, today, in any rigorous way. A possible approach to this problem is to turn to 
the phenomenological models of hadronization implemented in QCD-based parton-shower Monte Carlo 
(PSMC) codes. In PSMC the full final-state kinematical configuration is available at both the parton and 
hadron levels. Therefore, it is possible to "measure" D^\z) using eq. (16), both d<7# and da being 



known. In the present section, we carry out this program using the PSMC HERWIG [|28|]. HERWIG 
evolves quarks according to perturbative QCD down to small scales. The quarks are paired up at the 
end of the evolution into colour singlet clusters, which are then decayed to the physical hadrons using 
phenomenological distributions. The study of the heavy quark hadronisation process in HERWIG will 
allow us to test the universality assumption, and to measure the size of possible deviations. 

We should stress that, at this moment, our conclusions are only relevant for the hadronization 
model implemented in HERWIG; other PSMC's, which treat the hadronization process differently (for 
example, by adopting a string model), may well lead to different conclusions. 

In order to precisely define our procedure for extracting D^\z), we need to consider in more 
details the way in which HERWIG generates events. Regardless of the type of initial-state particles, we 
can distinguish the following steps. 

• Hard subprocess: at this stage, the PSMC generates the kinematics for the basic 2^2 hard reac- 
tion. We denote the momentum of the 6-quark (or antiquark) as p\ ard . 

• Parton shower: the partons resulting from the hard subprocess undergo successive branchings, 
until their virtuality is smaller than a fixed cutoff value. We denote the momentum of the 6-quark 
at the end of this phase as p^ s . 

• Gluon splitting and cluster formation: the gluons present at the end of the shower are decayed into 
light-quark pairs. Colour-singlet, two-body clusters are formed, according to colour parenthood 
and closeness in the phase-space. If there exist one or more cluster whose mass is too large 
(relative to a given threshold), part of the cluster rest energy is transformed into new qq pairs, and 
new clusters are defined. In this process, energy-momentum is redistributed among the cluster 
elements, and the momentum of the 6-quark can therefore be modified with respect to p^ s . The 
momentum of the b quark after completion of the clustering process will be denoted by pf sp . 

• Cluster decay and hadron formation: the clusters decay into observable hadrons, according to the 
flavour and to tabulated mass spectra. We therefore obtain ^-flavoured hadrons, whose momentum 
we denote as ps- 

The hard subprocess and parton shower stages are based on perturbative QCD. Thus, we identify the 
predictions given by HERWIG at the end of the parton shower with the cross section a that appears in 



eq. (16). On the other hand, the gluon splitting and cluster decay stages do not contain QCD informa- 
tion, as they are performed according to a phenomenological model. The g—>qq splitting and the decay 
kinematics are induced by simple phase-space considerations. We thus identify these stages as the long- 
distance, non perturbative part of the process, which gives rise to D^\z). We therefore determine the 
fragmentation function by comparing the results for p^ s and ps, defining, on an event-by-event basis, the 
following variables: 

_ PB ■ PT _ _Eb+Pb 1 pT_ 

\Pb \ E b + K I 

where p ps = p ps / |p^ s |- Our conclusions will apply to both z\ and z^, thus, we will collectively 
denote them by z. In hadronic collisions, the momenta and energies have to be substituted by transverse 



momenta and transverse energies respectively. Our results are summarized in table |11]; we considered 



e + e collisions at y/S = 91.2 GeV and pp collisions at \/S = 1.8 TeV. In the table, we present four of 



Table 1 1 : Normalized Mellin moments of the 6-quark non-perturbative fragmentation function. Results are given for the case 
of e + e~ collisions at vS = 91.2 GeV and for pp collisions at y/S = 1.8 TeV. All numbers have a statistical accuracy of 
±0.01. 





e + e 


PP 




Mi 


M2 


M3 




Mi 


M2 


M3 


M4 




pT 


< 5 GeV 


0.87 


0.78 


0.71 


0.66 


0.95 


0.94 


0.94 


0.96 


10 < 


pT 


< 15 GeV 


0.92 


0.85 


0.80 


0.76 


0.85 


0.74 


0.66 


0.60 


20 < 


->ps 

Pb 


< 25 GeV 


0.92 


0.86 


0.81 


0.76 


0.83 


0.71 


0.63 


0.57 


30 < 


—ps 

Pb 


< 35 GeV 


0.92 


0.85 


0.80 


0.75 


0.82 


0.70 


0.62 


0.56 



the (normalized) Mellin moments of the z distribution, defined as follows: 



J dzz n D%\z) J J dzD$\z). (18) 



Usually, < z < 1. In the present case, as we will see, we can also have z > 1; thus, in eq. (18) the 
range of integration coincide with the support of D^\z). The Mellin moments appearing in table [ll] 
have been evaluated by considering bins in \p^ s \ (in the case of hadronic collisions, the momentum is the 
transverse one). In e + e~ collisions larger (smaller) values of \p^ 3 \ correspond to less (more) energy lost 
to gluons. In hadronic collisions larger (smaller) values of \p b ps \ are more likely to correspond to larger 

(smaller) values of the hard process momentum before evolution. In either case, dependence of (z) 
on \plj S \ signals therefore a departure from universality. 

By inspection of the table, we see that (z) is scale-independent to a very good extent (the 
situation appears to be slightly better in the case of e + e~ collisions), except for the very low p^ s region; 
this is what we should expect, since in that region the factorization theorem on which eq. (|^) is based 
is bound to fail. On the other hand, there seems to emerge a clear difference between the fragmentation 
functions extracted from e + e~ and pp "data", the latter being substantially softer than the former. The 
first moment, which is the average value of the fragmentation variable, changes by about 10%. This 
variation can change the rate of predicted 6-hadrons in hadronic collisions by almost 50%. 

This suggests that transporting to hadronic collisions the non-perturbative fragmentation functions 
obtained by fitting e + e~ data may not be correct. Of course, a much more detailed investigation on the 
subject is required before reaching a firm conclusion; however, this simple exercise of ours shows that 
universality should not be taken for granted. 

We now concentrate on the separate role played in the fragmentation process by the purely per- 
turbative evolution and by the non-perturbative gluon-splitting phase, before the cluster formation and 
decay take place. We shall confine ourselves to the case of e + e~ collisions. The variables relevant to our 
study are the following: 

1. Energy fraction retained during the perturbative evolution: 



{Pb I _ {Pb I 



"P s /77 I ^>hard 



Vs \pt 



(19) 



where \/S is the e + e CM energy. 
2. Energy fraction retained during the gluon-splitting: 



->qlsp 
Pb 

Zglsp = \pPS\ ' 




3. Energy fraction left after the perturbative evolution and the gluon-splitting: 



Zp S X Zgl S p 



->qlsp 
Pb 



s 



->qlsp 
Pb 



Pb 



hard 



(21) 



The left panel in fig. 22 shows the three distributions for b quarks at = 91.2 GeV. The solid 
histogram represents the distribution of z ps . The distribution has the shape of a Gribov-Lipatov, with no 
indication of a Sudakov turn-over at large z ps . The dotted line is the distribution in z. A strong deforma- 
tion of the purely perturbative curve is clearly seen. The dashed line corresponds to the z g i sp distribution. 
This is part of what the MC treats as a non-perturbative component of the fragmentation process. The 
peak of the dashed histogram at z g [ sp = 1 corresponds to events where the cluster containing the heavy 
quark does not need to be further split, while the tail corresponds to events where the invariant mass of 
the heavy-quark cluster is too large, and additional light-quark pairs have to be produced by hand. Notice 
that almost as much energy is lost during this non-perturbative phase, as is lost during the perturbative 
evolution. 

For comparison, we show the same set of curves for the evolution of the charm quark (right panel 
of fig. 22). Notice that while the effect of the perturbative evolution is to soften the quark spectrum 
relative to the 6-quark case, the amount of energy lost due to gluon splitting is similar ((z glsp ) = 0.82, 
as opposed to {z h glsp ) = 0.85). This is bizarre, since one expects the non-perturbative part to scale with 
1/m^. The same result is found for the fragmentation of the s quark (left panel of fig. 23). Here {z s gl ) 



is 0.81. Again, a violation of the expected 1/m^, scaling. 

Things improve for the top quark, whose distributions for y/S = 2 TeV are shown on the right 



panel of fig. 23. The gluon-splitting part has only a minor impact on the overall spectrum of the top 
quark. 

We are a bit bothered by the dominant role played by the gluon-splitting phase. By comparison, 
the next step in the evolution, namely the cluster formation and decay, plays only a minor role, as will 
be shown next. We would have anticipated that the cluster formation and decay should be the place 
where most of the non-perturbative physics should show up. This suggests that the thresholds for the 
perturbative evolution in the MC shold be lowered, so that the impact of the non-perturbative gluon 
splitting phase is reduced, and purely perturbative Sudakov effects can manifest themselves. 

We now turn again to the non-perturbative part of the fragmentation function. The most striking 



feature, that cannot be inferred from the simple study of Mellin moments as done in table 1 1 , is the 



presence of a double peak in the high-z region (see the left panel of fig. |24|). A first peak (which we will 
call peak A) is seen at z values around 0.97. A second peak (peak B) is at z = 1.01 (we have a z-bin size 
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Fig. 23: Fragmentation functions for s (left) and t (right) quarks, produced in e + e collisions at y/s = 91.2 (left) and 2000 
(right) GeV. 
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Fig. 24: B-meson non-perturbative fragmentation function in e + e at ■s/S = 91.2 GeV, for 30 GeV< \f>b\ <35 GeV. 
CLDIR=1, CLSMR=0. See the text for details. 



of 0.02. We verified that the events contributing to the peak B do not have z = 1 + e, i.e., the peak is not 
due to a roundoff error). The latter peak is higher than the former. 

The origin of this double peak can be traced back to the following facts. First, the momentum 
of the emitted B meson is very strongly correlated with the momentum of the b quark which enters 
the cluster. Therefore, the z distribution closely reflects the mass spectrum of the light hadron emitted, 
together with the B meson, in the cluster decay. Second, the peak B is almost entirely due to events 
where the cluster decay into B + it: at this peak, z > 1 because the mass of the pion is lighter than the 
mass of the lightest quarks in HERWIG. More in detail, we have observed the following facts. 

• The structure of the double peak is strongly influenced by the value taken by the two input pa- 
rameters CLDIR and CLSMR. If the default is used (CLDIR=1, CLSMR=0), the double peak 



is observed (see the plot on the left of fig.[24|). On the other hand, by setting CLSMR^O, the z 
distributions display a single peak (broader that the previous ones) at about z = 0.97. For small 
CLSMR values and large b momenta, a second peak at z = 1.01 tends to re-appear, although 
smaller than observed before. 

The double peak disappears also if one chooses CLDIR=0, as in older HERWIG versions. In this 
case, the z distributions peak at about z = 0.9, this peak being much broader that those obtained 



with CLDIR=1, regardless of the value of CLSMR. 
• We then set CLDIR=1 and CLSMR=0. For any given B meson, we looked for the parent cluster 
C = {bcq}, and for the parent bottom quark, bp (the parent quark is defined as in the HERWIG 

routine HWCHAD). We observed what follows. 

- Plotting the z distributions for the events with p\, c / p\, p (i.e. events where the original 

cluster was split), we see a single peak, at the same z value as for the peak A (solid line, plot 
on the right of fig. 24). 

- The z distributions for events such that pb c = pb p display again a double peak. The two 
peaks are at the same z values as peaks A and B, the latter one being by far dominant. 

- Selecting only events with pt, c = p^ p , we found that the peak at the position of B corresponds 
to those clusters decaying into a B meson and a it, while the peak at the position of A is 
relevant for all the other two-body decays (dotted and dashed lines respectively, plot on the 
right of fig.||). 

Overall, notice also that the amount of energy retained after the gluon-splitting phase is of the same size 
as that retained at the end of the full hadronization process, indicating that cluster formation and decay 
have a minor impact on the total amount of energy lost during the non-perturbative part of the evolution. 

We were also able to reproduce the previous findings with a very simple model. Given a momen- 
tum for a quark b, we generate randomly the momentum for a light quark q, to be combined with b into 
a cluster, which eventually decays into a B meson and a particle of given mass rap. The momentum of 
the quark q is allowed to have a (small) transverse momentum with respect to the direction of the quark 
b. After evaluating the cluster mass, we performed the decay in the rest frame of the cluster, either in a 
isotropic manner (thus mimicking the choice CLDIR=0), or by letting the momentum of the meson B 
to be parallel to that of the quark b (which corresponds to CLDIR=1 and CLSMR=0). In the latter case, 
depending upon the value of mp, we got a peak for z < 1 (if mp > m q ) or z > 1 (if mp < m q ). 

In conclusion, the z distributions we find when using HERWIG seem not to contain a lot of dynam- 
ical information, the most important features being those implemented in the cluster-decay routine. If the 
decay is not smeared out (CLDIR=0), we get a structure which is very difficult to reconcile with the idea 
of fragmentation we have from QCD. After smearing, the distribution still has a z > 1 tail which will be 
extremely difficult to fit with a function vanishing for z—*l. This problem is related to the fact that the 
mass of the lightest quarks in the MC is 320 MeV, that is much larger than the pion mass. We performed 
a test by reducing the light quark masses to 20 MeV, and increasing the shower cutoff VQCUT in such 
a way as to maintain the default value of the effective infrared threshold. The double peak structure, as 
expected, disappeared. It remains to be seen, however, whether such a small value of the quark masses 
is, more generally, acceptable. 



4. A STUDY OF THE bb PRODUCTION MECHANISM IN PHYTIAf] 
4.1 Introduction 

In this section, we present a study on bb production performed within the CMS collaboration using the 
Monte Carlo package Pythia 5.75 as an event generator. In particular, we investigate the influence of 
the cut-off on the hard interaction transverse momentum P t on the production of bb events. 

In Monte Carlo programs, bb pairs in hadron collisions are produced by the mechanisms of gluon 
fusion, gluon splitting and flavour excitation. All of them give contributions of the same order to the total 
cross section, but they give rise to different kinematical configurations of the final state. 

There are two ways to generate bb events in Pythia: 

• Using a steering card MSEL=5, a gluon fusion mechanism (gg — ► bb) is mainly simulated. Each 
event contains at least one bb pair. 
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• Using MSEL=1, all QCD 2 — > 2 processes are simulated. In this case, all production mechanisms 
contribute to the bb production, but the probability to find a bb pair in the event is less than 1%. 

About one million events have been simulated in CMS with MSEL=1, in order to have a sample with 
all bb production mechanisms and default PYTHlAcut-off, not to introduce any bias in the kinematics. 
The selection efficiency of triggered events out of this sample is quite low. In order to have higher signal 



statistics, in some cases one can use kinematical cuts which are different from the PYTHlAdefault. [29] 



4.2 bb production 

Two samples have been prepared to investigate the influence of the P t cut on the production of bb events. 
Both of them have been generated using MSEL= 1 and contain events with only one bb pair. Only events 
with P t > 10 GeV have been selected in the samples. The first sample (SAMPLE A) has been generated 
with the default P t cut of 1 GeV and only events with P t > 10 GeV were selected. The second sample 
(SAMPLE B) has been generated with P t > 10 GeV. In both samples the following processes contribute: 

gg -> qq (22) 
gqi -> gq% (23) 

99 ^ 99 ■ (24) 



bb pair is produced by gluon splitting g — > bb in initial or final state shower evolution (processes (22) 



to (24)) or in the hard interaction (process (22)). For both samples A and B, we have computed the bb 



production cross section 

Jtot ~^bb -tot flG\ 

a bb ~ J^t a ' < 25 ) 

where is the number of bb events with P t > 10 GeV, N tot is the total number of generated events, 
and a tot is the total cross section (given by Pythia). We find that 

• for sample A, cr*5* = 150 fib. The gluon fusion contribution is about 20 fib, while the gluon 



splitting contributions are ~ 30 fib and ~ 100 fib for processes ( J23D and (|24|) respectively; 
for sample B, <r*2*=257 fib. Gluon fusion and gluon splitting contributions are at the same level as 
in sample A. In this case, however, there are also contributions from the processes bg — > bg and 
bq — ► bq of about 1 10 fib. In the following we will call these contributions flavour excitation. 



Figure 25 illustrates the difference in the bb production cross sections due to the additional contribution 
of the flavour excitation mechanism in sample B. The effect has the following explanation. When the 
default P t cut-off is used, PYTHIA generates processes in the low energy approximation, i.e. there are 
no heavy quarks inside the parton distribution. This approach changes if one uses a different Pt cut-off: 
the parton distributions in this case include also b and c quarks. As a consequence, samples A and B are 
different in two respects: values of the cross sections, and set of production mechanisms. The difference 
in the cross section is not very important, because the results are usually normalized to the total bb cross 
section of 500 fib. On the other hand, the different production mechanisms could be more dangerous, 
as they can lead to different kinematical distributions, and therefore affect the efficiencies of physical 
selection. 



4.3 Kinematics 

The main kinematical parameters which define the signature of bb event are the transverse momenta and 
pseudorapidities of the b quarks, and the angular distance Acft between their directions in the transverse 
plane. The first two parameters have similar distributions in both samples. The A<p distribution is shown 
in fig. 26 for the three different mechanisms. For what concerns gluon splitting, the distribution is slightly 
peaked at small Acj>. The angle between the two 6-quarks produced by the gluon-fusion mechanism has 
a peak at Acft ~ tt, as expected, since in the process gg — > bb the b-quarks are produced back-to-back in 




the transverse plane. The last distribution corresponds to the flavour excitation production mechanism, 
for which the back-to-back topology is preferred. We can conclude that the total Acf> distributions of 
sample A and sample B are slightly different. Some care should be taken about this, as it could affect the 
estimated efficiency of selection cuts. 

4.4 Pythia 6.125 

We have studied the same problem using the new 6.125 version PYTHIA. We have generated two new 
samples A and B with Pythia 6.125, and we have found the following results: 

• Sample A: the bb production cross section is cr*^*=220 fib. Gluon fusion contributes ~ 47 fib and 
total gluon splitting gives 173 fib. 

• Sample B: bb production cross section is cr*^' =465 fib. In this case, gluon fusion is ~ 51 fib and 
total gluon splitting is ~ 193 fib. The contribution from the flavour excitation is about 221 fib. 

Also the way Pythia 6.125 generates bb pairs depends on the P t cut-off. The difference with PYTHIA 
5.75 values, are due to different total cross section in the two versions. 
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Fig. 26: Angle between the two b quarks in the transverse plane: the upper one is gluon splitting, the middle is gluon fusion, 
the last is related to flavour excitation. In the plots the values are not normalized. 



4.5 Interpretation 

Many of the features of Pythia illustrated in this section are easily explained^. It turns out that Pythia 
treats differently processses with a low and high p™ m . The limit is related to the scale of multiple 
interactions, which is fixed to 2 GeV in the older versions, and was made energy dependent in PYTHIA 6, 
being 3.2 GeV at the LHC energy. When p^ m is above this scale, the hard process is selected according 
to conventional matrix elements. Below this scale, the hardest interaction is instead taken from the naive 
jet cross section multiplied by a "Sudakov style" form factor, that represents the probability that higher 
Pt interactions did not take place in the rest of the event. Since this procedure implies the computation 
of all parton-parton scattering processes, the choice was made to exclude from it the incoming b and 
c components, to save time in the computation. This feature is no longer considered useful in modern 
times, the computers being much faster. Thus, in Pythia 6.138, also the b and c processes will be 

4 T. Sjostrand and E. Norrbin, private communication. 



implemented in the low px mode. 

The difference in the total cross section in Pythia 5.7 and 6.1 have a physical origin, since 6.1 
uses newer parton distributions that, according to HERA data, are more singular in the small x region. 

The authors of Pythia recommend the following procedure for the generation of b events. Parton 
fusion and flavour excitation can be generated separately; the relevant massive matrix elements are used 
for parton fusion, and one can go to the limit pr^O with this process. Gluon splitting cannot be generated 
separately: all hard processes must be generated, excluding parton fusion and flavour excitation, and 
one should look for the heavy flavour. Multiple interactions are there switched off, in order to avoid a 
double-counting of the jet cross section. This is adequate for the study of the b production properties, but 
clearly does not fully represent the structure of the underlying event. In future PYTHIA versions, when 
flavour excitation is included in the minimum bias machinery with multiple interaction, this latter should 
offer an almost equivalent alternative, but still without the correct mass treatment of the parton fusion 
process near threshold. Other limitations still remain from complex problems related to the treatment of 
beam remnants; therefore, flavour excitation is only enabled for the hardest interaction in the multiple- 
interaction scenario. 

A sample of commented code is included below. By using different flags (MEKIND=0,1,2) three 
samples will be generated: parton fusion, flavour excitation and gluon splitting. 

INTEGER KFINTMP (-40 : 40) 
C... Multiple interactions switched off 

MSTP (81) =0 

PARP (81) =0 .DO 

PARP (82) =0 .DO 
C... Maximum virtuality in ISR is PARP ( 67 ) *Q* *2 
C... This is the default in Pythia 6.137 

PARP (67) =1 .DO 
C... Choose heavy quark (bottom=5, charm=4) 

MASSIVE=5 
C... Helper variable 

HQMASS=PMAS (MASSIVE, 1) 
C... Choose the kind of heavy quark production: 
C... MEKIND is a local variable set to 0, 1 or 2 

IF (MEKIND==0 ) THEN ! Massive matrix elements 
MSEL=MASSIVE 

ELSE IF (MEKIND==1) THEN ! Flavour excitation 
MSEL=1 

CKIN (3) =HQMASS 
CKIN (5) =CKIN (3) 
ELSE IF (MEKIND==2 ) THEN ! Gluon splitting (ISR, FSR) 
MSEL=1 

CKIN (3) =HQMASS 
CKIN (5) =CKIN (3) 
END IF 

C... More restrictive cuts can be put here. 
C... Example, 100 events in total. 

NEVENTS=100 
C *** EVENT LOOP *** 

IF (MEKIND==1) NEVENTS=NEVENTS / 2 
C. . . . Loop over incoming partons 

DO ISIDE=1,2 



IF (MEKIND/=1 .AND. ISIDE==1) THEN 

GOTO 100 
ELSE IF (MEKIND==1 ) THEN 
C... Only for flavour excitation: 
C... Make backup copy of KFIN array 
DO IKF=-40,40 

KFINTMP ( IKF ) =KFIN (ISIDE, IKF) 
C... Remove all incoming partons: 
KFIN (ISIDE, IKF) =0 
END DO 

C... Select only b/bbar as incoming partons: 
KFIN (ISIDE, MASSIVE) =1 
KFIN (ISIDE, -MASSIVE) =1 
END IF 

DO IEV=1, NEVENTS 
C. . . Generate an event 
CALL PYEVNT 

C... For gluon splitting, remove events with HQ in the hard interaction 
C... to avoid double counting: 
IF (MEKIND==2 ) THEN 
DO 1=5, 8 

IF (ABS (K (I, 2) ) ==MASSIVE) GOTO 50 
END DO 
END IF 
C . . . Analysis . . . 
50 END DO 

C... Print statistics 

CALL PYSTAT (1) 
C... Restore KFIN matrix: 

IF (MEKIND==1 . AND . ISIDE==1 ) THEN 
DO IKF=-40, 40 

KFIN (ISIDE, IKF) =KFINTMP (IKF) 
END DO 
END IF 
10 END DO 



5. ASYMMETRIES^ 
5.1 Introduction 

Sizeable leading particle asymmetries between e.g. D and D + have been observed in several fixed 
target experiments J30|]. It is of interest to investigate to what extent these phenomena translate to bottom 
production and higher energies. No previous experiment has observed asymmetries for bottom hadrons 
due to limited statistics or other experimental obstacles. Bottom asymmetries are in general expected to 
be smaller than for charm because of the larger bottom mass, but there is no reason why they should be 



absent. In the fixed target experiment HERA-B, bottom asymmetries could very well be large [31] even 
at central rapidities, but the conclusion of the present study is that asymmetries at the LHC are likely to 
be small. In the following we study possible asymmetries between B and B hadrons at the LHC within 



the Lund string fragmentation model [32] and the intrinsic heavy quark model [33] 
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In the string fragmentation model [|34]], the perturbatively produced heavy quarks are colour con- 
nected to the beam remnants. This gives rise to beam-drag effects where the heavy hadron can be 
produced at larger rapidities than the heavy quark. The extreme case in this direction is the collapse of 
a small string, containing a heavy quark and a light beam remnant valence quark of the proton, into a 
single hadron. This gives rise to flavour correlations which are observed as asymmetries. Thus, in the 
string model, there can be coalescence between a perturbatively produced bottom quark and a light quark 
in the beam remnant producing a leading bottom hadron. 

There is also the possibility to have coalescence between the light valence quarks and bottom 
quarks already present in the proton, because the wave function of the proton can fluctuate into Fock 
configurations containing a bb pair, such as | uudbb) . In these states, two or more gluons are attached 
to the bottom quarks, reducing the amplitude by 0(a 2 ) relative to parton fusion p5|]. The longest-lived 
fluctuations in states with invariant mass M have a lifetime of 0{2P\ & \ i /M 2 ) in the target rest frame, 
where Pi a b is the projectile momenta. Since the comoving bottom and valence quarks have the same 
rapidity in these states, the heavy quarks carry a large fraction of the projectile momentum and can thus 
readily combine to produce bottom hadrons with large longitudinal momenta. Such a mechanism can 
then dominate the hadroproduction rate at large xp. This is the underlying assumption of the intrinsic 



heavy quark model Q33|], in which the wave function fluctuations are initially far off shell. However, they 



materialize as heavy hadrons when light spectator quarks in the projectile Fock state interact with the 



target Q36Q. 

In both models the coalescence probability is largest at small relative rapidity and rather low trans- 
verse momentum where the invariant mass of the Qq system is small, enhancing the binding amplitude. 
One exception is at very large p±, where the collapse of a scattered valence quark with a b quark from 
the parton shower is also possible, giving a further (small) source of leading particle asymmetries in the 
string model. 

5.2 Lund String Fragmentation 

Before describing the Lund string fragmentation model, some words on the perturbative heavy quark 
production mechanisms included in the Monte Carlo event generator PYTHlA[f37|] used in this study is in 
order. We study pp events with one hard interaction because events with no hard interaction are not ex- 
pected to produce heavy flavours and events with more than one hard interaction — multiple interactions 
— are beyond the scope of this initial study and presumably would not influence the asymmetries. After 
the hard interaction is generated, parton showers are added, both to the initial (ISR) and final (FSR) state. 
The branchings in the shower are taken to be of lower virtualities than the hard interaction introducing 
a virtuality (or time) ordering in the event. This approach gives rise to several heavy quark production 
mechanisms, which we will call pair creation, flavour excitation and gluon splitting. The names may 
be somewhat misleading since all three classes create pairs at g^QQ vertices, but it is in line with the 
colloquial nomenclature. The three classes are characterized as follows. 

Pair creation The hard subprocess is one of the two LO parton fusion processes gg^QQ or qq^QQ. 
Parton showers do not modify the production cross sections, but only shift kinematics. For in- 
stance, in the LO description, the Q and Q have to emerge back-to-back in azimuth in order to 
conserve momentum, while the parton shower allows a net recoil to be taken by one or several 
further partons. 

Flavour excitation A heavy flavour from the parton distribution of one beam particle is put on mass 
shell by scattering against a parton of the other beam, i.e. Qq— >Qq or Qg— >Qg. When the Q is 
not a valence flavour, it must come from a branching g^QQ of the parton-distribution evolution. 
In most current sets of parton-distribution functions, heavy-flavour distributions are assumed to 
vanish for virtuality scales Q 2 < uiq. The hard scattering must therefore have a virtuality above 
rnq. When the initial-state shower is reconstructed backwards [38], the g^QQ branching will 
be encountered, provided that Qq, the lower cutoff of the shower, obeys Qq < m,Q. Effectively 



the processes therefore become at least gq^QQq or gg^QQg, with the possibility of further 
emissions. In principle, such final states could also be obtained in the above pair-creation case, 
but the requirement that the hard scattering must be more virtual than the showers avoids double 
counting. 

Gluon splitting A g^QQ branching occurs in the initial- or final-state shower but no heavy flavours are 
produced in the hard scattering. Here the dominant QQ source is gluons in the final-state showers 
since time-like gluons emitted in the initial state are restricted to a smaller maximum virtuality. 
Except at high energies, most initial state gluon splittings instead result in flavour excitation, al- 
ready covered above. An ambiguity of terminology exists with initial-state evolution chains where 
a gluon first branches to QQ and the Q later emits another gluon that enters the hard scattering. 
From an ideological point of view, this is flavour excitation, since it is related to the evolution 
of the heavy-flavour parton distribution. From a practical point of view, however, we choose to 
classify it as gluon splitting, since the hard scattering does not contain any heavy flavours. 
In summary, the three classes above are then characterized by having 2, 1 or 0, respectively, heavy 
flavours in the final state of the LO hard subprocess. Another way to proceed is to add next-to-leading 
order (NLO) perturbative processes, i.e the O(a^) corrections to the parton fusion [^] [Q]. However, with 
our currently available set of calculational tools, the NLO approach is not so well suited for exclusive 
Monte Carlo studies where hadronization is added to the partonic picture. 

Flavour excitation and gluon splitting give significant contributions to the total b cross section at 
LHC energies and thus must be considered when this is of interest, see the following. However, NLO 
calculations probably do a better job on the total b cross section itself (while, for the lighter c quark, 
production in parton showers is so large that the NLO cross sections are more questionable). The shapes 
of single heavy quark spectra are not altered as much as the correlations between Q and Q when extra 
production channels are added. Similar observations have been made when comparing NLO to LO 
calculations [||] Likewise, asymmetries between single heavy quarks are also not changed much by 
adding further production channels, so for simplicity we consider only the pair creation process here. 

After an event has been generated at the parton level we add fragmentation to obtain a hadronic 
final state. We use the Lund string fragmentation model. Its effects on charm production were described 
in [32]. Here we only summarize the main points. 

In the string model, confinement is implemented by spanning strings between the outgoing par- 
tons. These strings correspond to a Lorentz-invariant description of a linear confinement potential with 
string tension k k 1 GeV/fm. Each string piece has a colour charge at one end and its anticolour at the 
other. The double colour charge of the gluon corresponds to it being attached to two string pieces, while 
a quark is only attached to one. A diquark is considered as being in a colour antitriplet representation, 
and thus behaves (in this respect) like an antiquark. Then each string contains a colour triplet endpoint, 
a number (possibly zero) of intermediate gluons and a colour antitriplet end. An event will normally 
contain several separate strings, especially at high energies where g— >qq splittings occur frequently in 
the parton shower. 

The string topology can be derived from the colour flow of the hard process with some ambiguity 
arising from colour-suppressed terms. Consider e.g. the LO process gg— >bb where two distinct colour 
topologies are possible. Representing the proton remnant by a u quark and a ud diquark (alternatively 
d plus uu), one possibility is to have the three strings b-ud, b-u and u-ud, fig. 27, and the other is 
identical except the b is instead connected to the ud diquark of the other proton because the initial state 
is symmetric. 

Once the string topology has been determined, the Lund string fragmentation model [Q] can 
be applied to describe the nonperturbative hadronization. To first approximation, we assume that the 
hadronization of each colour singlet subsystem, i.e. string, can be considered separately from that of 
all the other subsystems. Presupposing that the fragmentation mechanism is universal, i.e. process- 
independent, the good description of e + e~ annihilation data should carry over. The main difference 
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Fig. 27: Example of a string configuration in a pp collision, (a) Graph of the process, with brackets denoting the final colour 
singlet subsystems, (b) Corresponding momentum space picture, with dashed lines denoting the strings. 



between e + e and hadron-hadron events is that the latter contain beam remnants which are colour- 
connected with the hard-scattering partons. 

Depending on the invariant mass of a string, practical considerations lead us to distinguish the 
following three hadronization prescriptions: 

Normal string fragmentation In the ideal situation, each string has a large invariant mass. Then the 
standard iterative fragmentation scheme, for which the assumption of a continuum of phase-space 
states is essential, works well. The average multiplicity of hadrons produced from a string in- 
creases linearly with the string 'length', which means logarithmically with the string mass. In 
practice, this approach can be used for all strings above some cutoff mass of a few GeV. 

Cluster decay If a string is produced with a small invariant mass, perhaps only a single two-body final 
state is kinematically accessible. In this case the standard iterative Lund scheme is not applicable. 
We call such a low-mass string a cluster and consider its decay separately. When kinematically 
possible, a Q-q cluster will decay into one heavy and one light hadron by the production of a 
light qq pair in the colour force field between the two cluster endpoints with the new quark flavour 
selected according to the same rules as in normal string fragmentation. The q cluster end or the 
new qq pair may also denote a diquark. In the latest version of Pythia, anisotropic decay of a 
cluster has been introduced, where the mass dependence of the anisotropy has been matched to 
string fragmentation. 

Cluster collapse This is the extreme case of cluster decay, where the string mass is so small that the 
cluster cannot decay into two hadrons. It is then assumed to collapse directly into a single hadron 
which inherits the flavour contents of the string endpoints. The original continuum of string/cluster 
masses is replaced by a discrete set of hadron masses, mainly B and B* (or the corresponding 
baryon states). This mechanism plays a special role since it allows flavour asymmetries favouring 
hadron species that can inherit some of the beam-remnant flavour contents. Energy and momentum 
is not conserved in the collapse so that some energy-momentum has to be taken from, or transferred 
to, the rest of the event. In the new version, a scheme has been introduced where energy and 
momentum are shuffled locally in an event. 

We assume that the nonperturbative hadronization process does not change the perturbatively cal- 



culated total rate of bottom production. By local duality arguments Q39Q, we further presume that the rate 



of cluster collapse can be obtained from the calculated rate of low-mass strings. In the process e + e - ^cc 
local duality suggests that the sum of the 3/ip and ip' cross sections approximately equal the perturbative 
cc production cross section in the mass interval below the DD-threshold. Similar arguments have also 



been proposed for r decay to hadrons [40] and shown to be accurate. In the current case, the presence 



of other strings in the event also allows soft-gluon exchanges to modify partem momenta as required 



to obtain the correct hadron masses. Traditional factorization of short- and long-distance physics would 
then also preserve the total bottom cross section. Local duality and factorization, however, do not specify 
how to conserve the overall energy and momentum of an event when a continuum of bd masses is to be 
replaced by a discrete B°. In practice, however, the different possible hadronization mechanisms do not 
affect asymmetries much. The fraction of the string-mass distribution below the two particle threshold 
effectively determines the total rate of cluster collapse and therefore the asymmetry. 

The cluster collapse rate depends on several model parameters. The most important ones are listed 
here with the Pythia parameter values that we have used. The Pythia parameters are included in the 
new default parameter set in Pythia 6.135 and later versions. 

• Quark masses The quark masses affect the threshold of the string-mass distribution. Changing 
the quark mass shifts the string-mass threshold relative to the fixed mass of the lightest two-body 
hadronic final state of the cluster. Smaller quark masses imply larger below-threshold production 
and an increased asymmetry. The new default masses are PMAS ( 1 ) = m u = PMAS ( 2 ) = = 
0.33D0, PMAS ( 3 ) = m s = 0.5D0, PMAS ( 4 ) = m c = 1.5D0 and PMAS ( 5 ) = m b = 4.8D0. 

• Width of the primordial k± distribution. If the incoming partons are given small p± kicks in the 
initial state, asymmetries can appear at larger pj_ since the beam remnants are given compensating 
p± kicks, thus allowing collapses at larger p±. The new parameters are PARP ( 91 ) =1 . DO and 
PARP (93)=5.D0. 

• Beam remnant distribution functions (BRDF). When a gluon is picked out of the proton, the rest 
of the proton forms a beam remnant consisting, to first approximation, of a quark and a diquark. 
How the remaining energy and momentum should be split between these two is not known from 
first principles. We therefore use different parameterizations of the splitting function and check 
the resulting variations. We find significant differences only at large rapidities where an uneven 
energy-momentum splitting tend to shift bottom quarks connected to a beam remnant diquark more 
in the direction of the beam remnant, hence giving rise to asymmetries at very large rapidities. We 
use an intermediate scenario in this study, given by MSTP ( 92 ) =3. 

• Threshold behaviour between cluster decay and collapse. Consider a bd cluster with an in- 
variant mass at, or slightly above, the two particle threshold. Should this cluster decay to two 
hadrons or collapse into one? In one extreme point of view, a Bn pair should always be formed 
when above this threshold, and never a single B. In another extreme, the two-body fraction would 
gradually increase at a succession of thresholds: Bit, B*tt, Bp, B*p, etc., where the relative prob- 
ability for each channel is given by the standard flavour and spin mixture in string fragmentation. 
In our current default model, we have chosen to steer a middle course by allowing two attempts 
(MST J ( 1 7 ) =2) to find a possible pair of hadrons. Thus a fraction of events may collapse to a 
single resonance also above the Bit threshold, but Btt is effectively weighted up. If a large number 
of attempts had been allowed (this can be varied using the free parameter MST J ( 1 7 ) ), collapse 
would only become possible for cluster masses below the Bit threshold. 

The colour connection between the produced heavy quarks and the beam remnants in the string 
model gives rise to an effect called beam remnant drag. In an independent fragmentation scenario the 
light cone energy momentum of the quark is simply scaled by some factor picked from a fragmentation 
function. Thus, on average the rapidity is conserved in the fragmentation process. This is not necessarily 
so in string fragmentation, where both string ends contribute to the four-momentum of the produced 
heavy hadron. If the other end of the string is a beam remnant, the hadron will be shifted in rapidity in 
the direction of the beam remnant resulting in an increase in \y\. This beam-drag is shown qualitatively 
in fig. 28, where the rapidity shift is shown as a function of rapidity and transverse momentum. This 
shift is not directly accessible experimentally, only indirectly as a discrepancy between the shape of 
perturbatively calculated quark distributions and the data. 
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Fig. 28: (a) Average rapidity shift Ay = (ye — J/b) as a function of y for some different p± cuts, (b) Average rapidity shift 
(| Ay\) in the direction of "the other end of the string" that the bottom quark is connected to, i.e. ignoring the sign of the shift. 



5.3 Intrinsic Heavy Quarks 

The wavefunction of a hadron in QCD can be represented as a superposition of Fock state fluctuations, 
e.g. \ny), \nyg), |ravQQ), ■ ■ ■ components where ray = uud for a proton. When the projectile scatters 
in the target, the coherence of the Fock components is broken and the fluctuations can hadronize either 
by uncorrected fragmentation as for leading twist production or coalescence with spectator quarks in the 
wavefunction [ ]33| ] p6[]. The intrinsic heavy quark Fock components are generated by virtual interactions 
such as gg — > QQ where the gluons couple to two or more projectile valence quarks. Intrinsic QQ 
Fock states are dominated by configurations with equal rapidity constituents so that, unlike sea quarks 
generated from a single parton, the intrinsic heavy quarks carry a large fraction of the parent momentum 

The frame-independent probability distribution of an n-particle bb Fock state is 

dPl = N Si) (26) 

dxi ■ ■ ■ dx n n (m 2 - ££=i("»i/zi)) 2 ' 

where fhf = k\ i + mf is the effective transverse mass of the i th particle and X{ is the light-cone mo- 
mentum fraction. The probability, Pfc, is normalized by N n and n = 5 for baryon production from the 
Iriybb) configuration. The delta function conserves longitudinal momentum. The dominant Fock con- 
figurations are closest to the light-cone energy shell and therefore the invariant mass, M 2 = ^ fhf/xi, 
is minimized. Assuming {k 2 Li ) is proportional to the square of the constituent quark mass, we choose 
m q = 0.45 GeV, m s = 0.71 GeV, and fh h = 5 GeV [|l|] @]. 

The xp distribution for a single bottom hadron produced from an n-particle intrinsic bottom state 
can be related to P£ and the inelastic pp cross section by 



^(PP) _ f^ffin Jf!_„4 

dx F ~ dx F ™4ml sK1 ""' 



< p -^af(M bE ) . (27) 



The probability distribution is the sum of all contributions from the |nvbb) and the |nybbqq) configura- 
tions with q = u, d, and s and includes uncorrected fragmentation and coalescence, as described below, 
when appropriate The factor of u 2 /4m^ arises from the soft interaction which breaks the coherence 
of the Fock state. We take u 2 ~ 0.1 GeV 2 ||44|]. The intrinsic charm probability, P^ c = 0.31%, was deter- 
mined from analyses of the EMC charm structure function data p5|]. The intrinsic bottom probability is 
scaled from the intrinsic charm probability by the square of the transverse masses, = P\ c (fh c / fh^) 2 . 



The intrinsic bottom cross section is reduced relative to the intrinsic charm cross section by a factor of 
a^{M h -^) / (M c5 ) @]. Taking these factors into account, we obtain af h (pN) « 7 nb at 14 TeV. 

There are two ways of producing bottom hadrons from intrinsic bb states. The first is by un- 
corrected fragmentation. If we assume that the b quark fragments into a B meson, the B distribution 
is 

= dz\ \ dxi- ^ '- 6{xb - zx h ) , (28) 

dx-Q J "7 dx\ . . . dx n z 

These distributions are assumed for all intrinsic bottom production by uncorrected fragmentation with 
D H / b (z) = 5(z — 1). At low p±, this approximation should not be too bad, as seen in fixed target 
production [|42|]. 

If the projectile has the corresponding valence quarks, the bottom quark can also hadronize by 
coalescence with the valence spectators. The coalescence distributions are specific for the individual 
bottom hadrons. It is reasonable to assume that the intrinsic bottom Fock states are fragile and can easily 
materialize into bottom hadrons in high-energy, low momentum transfer reactions through coalescence. 
The coalescence contribution to bottom hadron production is 



1=1 

where the coalescence function is simply a delta function combining the momentum fractions of the 
quarks in the Fock state configuration that make up the valence quarks of the final-state hadron. 

Not all bottom hadrons can be produced from the minimal intrinsic bottom Fock state configu- 
ration, |nybb). However, coalescence can also occur within higher fluctuations of the intrinsic bottom 
Fock state. For example, in the proton, the B and H£ can be produced by coalescence from |nvbbuu) 
and |nybbss) configurations. These higher Fock state probabilities can be obtained using earlier re- 
sults on tpip pair production [47] [S3]. If all the measured tpip pairs [ |49| ] arise from |nvcccc) config- 



urations, Pi cc « 4.4% Pi c [48] [|50|]. It was found that the probability of a |nvccqq) state was then 
Pcq = (^c/"iq) 2 Picc [ 47 ]. If we then assume Pb q = (m- c /^b) 2 Picq> we find that 

^ibq » ^ M P «* > (3°) 

\m h J \m q J 

leading to P^ u = P;bd ~ 70.4% P^ and P ic , s ~ 28.5% P^. To go to still higher configurations, one 
can make similar assumptions. However, as more partons are included in the Fock state, the coalescence 
distributions soften and approach the fragmentation distributions, eventually producing bottom hadrons 
with less momentum than uncorrected fragmentation from the minimal bb state if a sufficient number 
of qq pairs are included. There is then no longer any advantage to introducing more light quark pairs 
into the configuration — the relative probability will decrease while the potential gain in momentum is 
not significant. Therefore, we consider production by fragmentation and coalescence from the minimal 
state and the next higher states with uu, dd and ss pairs. 



The probability distributions entering Eq. (27) for B and B are 



dP B0 _ i/i dPl F idPl c \ i/i dPl F ldPZS 



\ Ui dPZ idp^\ 



dxF 2 \10 dxF 4 dxp ) 2 \10 dxF 5 dxF J 

, 1 M . . 1 + 1 ^Sl) (3l) 

2\10 dx F 5 dx F J 2\10 dx F 5 dx F J 

<ipv i dp$? , i dPZ , i ( i dPZ , idPi%\ , i dpg 



_|_ _ j ipq | ipq 1 _|_ ids (32) 

dxF 10 dxp 10 dxF 2 I 10 dx-p 8 dxp J 10 dxp 



See Ref. [43] for more details and the probability distributions of other bottom hadrons. 



5.4 Model predictions 



o 



In this section we present some results from both models. Figure |29| shows the asymmetry between B 
and B° as a function of y for several p± cuts in the string model. The asymmetry is essentially zero 
for central rapidities and increases slowly with rapidity. When the kinematical limit is approached, the 
asymmetry changes sign for small p± because of the drag-effect since b-quarks are often connected to 
diquarks from the proton beam remnant, fig. 27, thus producing B hadrons which are shifted more in 
rapidity than B°. Cluster collapse, on the other hand, tend to enhance the production of leading particles 
(in this case B°) so the two mechanisms give rise to asymmetries with different signs. Collapse is the 
main effect at small rapidities while eventually at very large y, the drag effect dominates. 
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Fig. 29: The asymmetry, A — "^go^^g" j , as a function of rapidity for different p± cuts: (a) p± < 5, 10 GeV and (b) 
p± > 5, 10 GeV using parameter set 1 as described in the text. 



In Table 12 we study the parameter dependence of the asymmetry by looking at the integrated 



asymmetry for different kinematical regions using three different parameter sets: 



• Set 1 is the new default as presented in section p.2 . 

• Set 2 The same as Set 1 except it uses simple counting rules in the beam remnant splitting, i.e. 
each quark get on average one third of the beam remnant energy-momentum. 

• Set 3 The old parameter set, before fitting to fixed-target data, is included as a reference. This 
set is characterized by current algebra masses, lower intrinsic k±, and an uneven sharing of beam 
remnant energy-momentum. 

We see that in the central region the asymmetry is generally very small whereas for forward (but 
not extremely forward) rapidities and moderate pj_ the asymmetry is around 1-2%. In the very for- 
ward region at small p±, drag asymmetry dominates which can be seen from the change in sign of the 
asymmetry. The asymmetry is fairly stable under moderate variations in the parameters even though the 
difference between the old and new parameter sets (Set 1 and 3) are large in the central region. Set 1 
typically gives rise to smaller asymmetries. 

The cross sections for all intrinsic bottom hadrons are given as a function of xp in fig. |3(]. The 
bottom baryon distributions are shown in fig. 3C(a). The A° distributions are the largest and most 



forward peaked of all the distributions. The E^" is the smallest and the softest, similar to that of the 
bottom-strange mesons and baryons shown in fig. |3^(b). The different coalescence probabilities as- 
sumed for hadrons from the |uudbbss) configuration have little real effect on the shape of the cross 



section, dominated by independent fragmentation. Of the B mesons shown in fig. 30 [c), the B + and B 



)0 



cross sections are the largest since both can be produced from the 5 particle configuration. The B and 
B° distributions are virtually identical. We note that the xf distributions of other bottom hadrons not 



Table 12: Parameter dependence of the asymmetry in the string model. The statistical error in the last digit is shown in 
parenthesis (95% confidence). 



Parameters 


\y\ < 2.5, p ± > 5GeV 


3 < |y| < 5, p± > 5 GeV 


\y\ >3, P± <5 GeV 


Set 1 


0.003(1) 


0.015(2) 


-0.008(1) 


Set 2 


-0.000(2) 


0.009(3) 


-0.005(2) 


Set 3 


0.013(2) 


0.020(3) 


-0.018(2) 



included in the figure would be similar to the bottom-strange hadrons since they would be produced by 
fragmentation only. 
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Fig. 30: Predictions for bottom hadron production are given for pp collisions at 14 TeV. The bottom baryon distributions are 
given in (a) for A° = E° (dot-dashed), SjJ" (dashed), and E^" (solid). The bottom-strange distributions are shown in (b) for H° 
(solid), (dashed), B|? (dot-dashed), and (dotted). In (c), the B meson distributions are given: B + (solid), B~ (dashed), 
B° (dot-dashed), and B° (dotted). The B~ and B° distributions are virtually identical. 

The xf distribution for final-state hadron H is the sum of the leading-twist fusion and intrinsic 
bottom components, 



da hN 



da* + dag 



dxp dxp dxF 



(33) 



The intrinsic bottom cross sections from Section 5.3 are combined with a leading twist calculation using 
independent fragmentation where drag effects are not included. The leading twist results have been 



smoothed and extrapolated to large xp to facilitate a comparison with the intrinsic bottom calculation. 
The resulting total B° and B° distributions are shown in fig. along with the corresponding asymmetry. 
Note that since the intrinsic heavy quark p± distributions are more steeply falling than the leading twist, 
we only consider < 5 GeV. The distributions are drawn to emphasize the high xp region where the 
distributions differ. The asymmetry is ~ 0.1 at xp ~ 0.25, corresponding to y ~ 6.5. Therefore, intrinsic 
bottom should not be a significant source of asymmeUies. 
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Fig. 31: (a) Leading-twist predictions for B° (solid) and B° (dashed) using independent fragmentation. Model predictions for 
B° (dot-dashed) and B° (dotted) distributions from Eq. (^3|). (b) The asymmetry between B° and B°, the dot-dashed and dotted 
curves in (a), is also given. 



5.5 Summary 

To summarize, we have studied possible production asymmetries between b and b hadrons, especially 
B° and B , as predicted by the Lund string fragmentation model and the intrinsic heavy quark model. 
We find negligible asymmetries for central rapidities and large p± (in general, less than 1%). For some 
especially favoured kinematical ranges such as y > 3 and 5 < p± < 10 GeV the collapse asymmetry 
could be as high as 1-2%. Intrinsic bottom becomes important only for xp > 0.25 and p± < 5 GeV, 
corresponding to y > 6.5. 



6. QUARKONIUM PRODUCTION^ 

The production of charmonium and bottomonium states at high-energy colliders has been the subject of 
considerable interest during the past few years. New experimental results from pp, ep and e + e~ colliders 
have become available, some of which revealed dramatic shortcomings of earlier quarkonium production 
models. In theory, progress has been made on the factorization between the short distance physics of 
heavy-quark creation and the long-distance physics of bound state formation. The colour-singlet model 
[ p| ] p2[ ] has been superseded by a consistent and rigorous framework, based on the use of non-relativistic 
QCD (NRQCD) [53], an effective field theory that includes the so-called colour-octet mechanisms. On 
the other hand, the colour evaporation model 



[55] [56] of the early days of quarkonium physics has 
been revived [57] [ j58| ] [59] [pQ]. However, despite the recent theoretical and experimental developments 
the range of applicability of the different approaches is still subject to debate, as is the quantitative 
verification of factorization. Because the quarkonium mass is still not very large with respect to the QCD 



scale, in particular for the charmonium system, non-factorizable corrections [61] 0621] [ 63 ] may not be 



" Thanks to J. Klay at UC Davis for extending the curves to large if. 
6 Section coordinators: M. Kramer, F. Maltoni, M.A. Sanchis-Lozano 



suppressed enough, if the quarkonium is not part of an isolated jet, and the expansions in NRQCD may 
not converge very well. In this situation a global analysis of various processes is mandatory in order 
to assess the importance of different quarkonium production mechanisms, as well as the limitations 
of a particular theoretical framework (for reviews on different quarkonium production processes see 
e.g. [|4|] [65] [|66[].) By the time the LHC starts operating, new experimental data from the Tevatron 



and HERA as well as theoretical progress, e.g. in the calculation of higher-order corrections, will have 
significantly improved the present picture and will allow more precise predictions than what is possible 
at present. In the following, we will therefore focus on the general phenomenological implications 
of the NRQCD approach for quarkonium production at the LHC, rather than aiming at a detailed and 
comprehensive numerical analysis. Based on the information provided by the present Tevatron data we 
will derive predictions for observables crucial for future LHC analyses, such as differential cross sections 
and quarkonium polarization. 

In the NRQCD approach, the cross section for producing a quarkonium state H at a hadron collider 
can be expressed as a sum of terms, each of which factors into a short-distance coefficient and a long- 
distance matrix element: 

da(pp/pp^H + X) = Y^ da(pp/pp^QQ [n] + x) (O h [n]> , (34) 

n 

where n denotes the colour, spin and angular momentum state of an intermediate QQ pair. The short- 
distance cross section da can be calculated perturbatively in the strong coupling a s . The NRQCD matrix 
elements (O [n]) (see [53] for their definition) are related to the non-perturbative transition probabilities 



from the QQ state n into the quarkonium H . They scale with a definite power of the intrinsic heavy-quark 
velocity v [67]. (v 2 ~ 0.3 for charmonium and v 2 ~ 0.1 for bottomonium.) The general expression (34) 



is thus a double expansion in powers of a s and v. 

The NRQCD formalism implies that so-called colour-octet processes associated with higher Fock 
state components of the quarkonium wave function must contribute to the cross section. Heavy quark 
pairs that are produced at short distances in a colour-octet state can evolve into a physical quarkonium 
through radiation of soft gluons at late times in the production process, when the quark pair has already 
expanded to the quarkonium size. Such a possibility is ignored in the colour-singlet model, where only 
those heavy quark pairs that are produced in the dominant Fock state (i.e. in a colour-singlet state and 
with the spin and angular momentum quantum numbers of the meson) are assumed to form a physical 
quarkonium. The most profound theoretical evidence that the colour-singlet model is incomplete comes 
from the presence of infrared divergences in the production cross sections and decay rates of P-wave 
states. Within the NRQCD approach, this problem finds its natural solution since the infrared singu- 



larities are factored into a colour-octet operator matrix element [68]. While colour-octet contributions 
are needed for a consistent description of P-wave quarkonia, they are phenomenologically even more 
important for 5- wave states like J ftp or T. According to the velocity scaling rules, colour-octet ma- 
trix elements for the production of 5-wave quarkonia are suppressed by a factor v 4 compared to the 
leading colour-singlet contributions. However, as discussed in some detail below, colour-octet processes 
can become significant if the short-distance cross section for producing QQ in a colour-octet state is 
enhanced. 

The production of 5-wave charmonium in pp collisions at the Tevatron has attracted considerable 
attention and has stimulated much of the recent theoretical development in quarkonium physics. The 
CDF collaboration has measured cross sections for the production of J/ip and ip(2S) states not coming 



from B or radiative \ decays, for a wide range of transverse momenta 5GeV <, pt(ip) <> 20GeV [69] 
[|7C|]. Surprisingly, the experimental cross sections were found to be orders of magnitudes larger than the 
theoretical expectation based on the leading-order colour-singlet model j f7"T| ] y72[]. This result is partic- 
ularly striking because the data extends out to large transverse momenta where the theoretical analysis 
is rather clean. The shortcoming of the colour-singlet model can be understood by examining a typical 
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Fig. 32: Generic diagrams for J/ip production in hadron-hadron collisions via colour-singlet and colour-octet channels. 



Feynman diagram contributing to the leading-order parton cross section, fig. ^2|(a). At large transverse 
momentum, the two internal quark propagators are off-shell by ~ pi so that the parton differential cross 
section scales like da /dpi ~ 1/pf. as indicated in the figure. On the other hand, when p t 3> 2m c the 
quarkonium mass can be considered small and the inclusive charmonium cross section is expected to 
scale like any other single-particle inclusive cross section ~ l/pf. The dominant production mechanism 



for charmonium at sufficiently large p t must thus be via fragmentation [73], the production of a parton 



with large p t which subsequently decays into charmonium and other partons. A typical fragmentation 
contribution to colour-singlet J /ip production is shown in fig. |32|(b). While the fragmentation contribu- 
tions are of higher order in a s compared to the fusion process fig. |32|(a), they are enhanced by a power 
pf/(2m c ) A at large p t and can thus overtake the fusion contribution at p t 3> 2m c . When colour-singlet 
fragmentation is included, the p t dependence of the theoretical prediction is in agreement with the Teva- 
tron data but the normalization is still underestimated by about an order of magnitude [[74]] [75] 



indicating that an additional fragmentation contribution is still missing. It is now generally believed that 
gluon fragmentation into colour-octet 3 Si charm quark pairs [ [77| ] J78|], as shown in fig. f32|(c), is the 
dominant source of J/ijj and tp(2S) at large p t at the Tevafron. The probability of forming a J/ijj particle 
from a pointlike cc pair in a colour-octet 3 Si state is given by the NRQCD matrix element (O J W [ 3 s{ 8) ]} 
which is suppressed by v A relative to the non-perturbative factor of the leading colour-singlet term. How- 
ever, this suppression is overcompensated for by the gain in two powers of a s /ir in the short-distance 
cross section for producing colour-octet 3 Si charm quark pairs as compared to colour-singlet fragmen- 
tation. At 0(v ) in the velocity expansion, two additional colour-octet channels have to be included, 



fig. 32(d), which do not have a fragmentation interpretation at order a 3 but which become significant at 



moderate p t ~ 2m c [79] [80]. The importance of the and 3 Pj 8 ' contributions cannot be estimated 



3p(8) 



from naive power counting in a s and v alone, but rather follows from the dominance of t-channel gluon 
exchange, forbidden in the leading-order colour-singlet cross section. 

The different contributions to the J/ ip transverse momentum distribution are compared to the CDF 
data [70] in fig. [53]. As mentioned above, the colour-singlet model at lowest order in a s fails dramati- 
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Fig. 33: Colour-singlet and colour-octet contributions to direct J/tp production inpp^J/ip+X at the Tevatron (y/s = 1.8 TeV, 
pseudorapidity cut |r/| < 0.6)) together with experimental data from CDF [|7C|]. Parameters: CTEQ5L parton distribu- 
tion functions jj8l]|; factorization and renormalization scale /i = \/ pf + 4mj?; m c — 1.5 GeV. The leading logarithms 
(a s lnp? /(2m c ) 2 ) n have been summed by solving the Altarelli-Parisi evolution equations for the gluon fragmentation func- 
tion. NRQCD matrix elements as specified in Table [l3|. 



Table 13: NRQCD matrix elements for charmonium production. The colour-singlet matrix elements are taken from the po- 
tential model calculation of [ [s^t j83|]. The colour-octet matrix elements have been extracted from the CDF data [|7c|], where 
MkC S o\ 3p o a) ) = (Os^So]) + k(Og [ 3 P ])/m 2 . The errors quoted are statistical only. Parameters: CTEQ5L parton 
distribution functions renormalization and factorization scale /j, = (p 2 +4m 2 l ) 1 / 2 an d m c = 1.5 GeV. The Altarelli-Parisi 
evolution has been included for the 3 S^ 8 ' fragmentation contribution. See [B] for further details. 
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cally when confronted with the experimental results. When colour-singlet fragmentation is included, the 
prediction increases by more than an order of magnitude at large p t , but it still falls below the data by a 
factor of ~ 30. The CDF results on charmonium production can be explained by including the leading 
colour-octet contributions and adjusting the unknown non-perturbative parameters to fit the data. Nu- 
merically one finds the non-perturbative matrix elements to be of 0(1O -2 GeV 3 ), see Table |l3|, perfectly 
consistent with the v 4 suppression expected from the velocity scaling rules. Similar conclusions can be 
drawn for ip(2S) production at the Tevatron. 

The analysis of the CDF data alone, although very encouraging, does not strictly prove the phe- 
nomenological relevance of colour-octet contributions because free parameters have to be introduced to 
fit the data. However, if factorization holds the non-perturbative matrix elements, Table 13, are universal 
and can be used to make predictions for various processes and observables. Besides a global analysis of 
different reactions, the measurement of quarkonium cross sections at the LHC will be crucial to assess 
the importance of the individual production mechanisms and to test factorization. In fig. 34 we have 
collected the cross section predictions for direct J/tp and ip(2S) production as well as the production 
of J/ip from radiative x decays at the LHC. The theoretical curves include the statistical errors in the 
extraction of the NRQCD matrix elements [Table 13]. There are, however, additional theoretical uncer- 
tainties which might affect the prediction, but which have not yet been fully quantified. In particular 
the determination of the (Og [ 1 5o]) and (Of [ 3 Po]} matrix elements (ip denoting J/ 1/3 or ip(2S)) from 
the Tevatron data is very sensitive to effects that modify the shape of the charmonium p t distribution at 
relatively small p t <^ 8 GeV. Those effects include the small-x behaviour of the gluon distribution [84], 
the evolution of the strong coupling [p4f|, as well as systematic effects inherent in NRQCD, such as the 
inaccurate treatment of the energy conservation in the hadronization of the colour-octet cc pairs [85]. 
Moreover, higher-order QCD corrections are expected to play an important role, as discussed in more 
detail below. The cross sections collected in fig. |4] should thus not be viewed as firm NRQCD predic- 
tions but will be refined as more experimental and theoretical information on charmonium production 
becomes available over the next few years. 

The inclusion of higher-order QCD corrections is required to reduce the theoretical uncertainty 
and to allow a more precise prediction of the LHC cross sections. Next-to-leading order (NLO) calcula- 
tions for quarkonium production at hadron colliders are presently available only for total cross sections 
[ |86| ] [87]. Significant higher-order corrections to differential distributions are expected from the strong 
renormalization and factorization scale dependence of the leading-order results [p5Q. Moreover, the NLO 
colour-singlet cross section includes processes like g + g—>QQ[ s S^] + g + g which are dominated by 
t-channel gluon exchange and scale as ~ a\(2mQ) 2 jp\. At p t S> 2toq their contribution is enhanced 
with respect to the the leading-order cross section, fig. |32](a), which scales as ~ a 3 (2mg) 4 /pf . This is 
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Fig. 34: Cross sections for J/ip and ip(2S) production in pp^tp + X at the LHC {-^/s — 14 TeV, pseudorapidity cut \r)\ < 2.5). 
Parameter specifications as in fig. The leading logarithms (a s lnpf / (2m c ) 2 )" have been summed by solving the Altarelli- 
Parisi evolution equations for the gluon fragmentation function. The error bands include the statistical errors in the extraction 
of the NRQCD matrix elements [Table |l|] only. 



born out by preliminary studies [ |88| ] which include part of the NLO hadroproduction cross section and 
by the complete calculation of NLO corrections to the related process of quarkonium photoproduction 
[j39|]. The NLO colour-singlet cross section may be comparable in size to the colour-octet 1 Sq and 3 Pj 
processes, which scale as ~ a^v 4 (2mg) 2 / 'pf (see fig. 32 [d)), and affect the determination of the cor- 
responding NRQCD matrix elements from the Tevatron data. A full NLO analysis is however needed 
before quantitative conclusions can be drawn. 

Another source of potentially large higher-order corrections is the multiple emission of soft or 
almost collinear gluons from the initial state partons. These corrections, as well as effects related to in- 
trinsic transverse momentum, are expected to modify the shape of the transverse momentum distribution 
predominantly at relatively low values of p t <, 2mQ. Initial state radiation can be partially summed in 



perturbation theory [gOp, but so far only total cross sections have been considered in the literature 
An estimate of the effect on the transverse momentum distribution should be provided by phenomeno- 
logical models where a Gaussian /c t -smearing is added to the initial state partons. The result of these 
calculations not only depends on the average {kt}, which enters as a free parameter, but also on the de- 
tails of how the smearing is implemented. Moreover, a lower cut-off has to be provided which regulates 
the divergences at p t = 0. Using the NLO calculation for the total cross section [|87|], one can obtain the 
rough estimate that perturbative Sudakov effects should be confined below p t ~ 1 — 2 GeV for both char- 
monium and bottomonium production at Tevatron energies. Qualitatively, the inclusion of /c t -smearing 
leads to an enhancement of the short distance cross section at small p t , which results in smaller values 



for the fits of the {C$1 1 S ]) and (O$[ 3 P ]) NRQCD matrix elements [||] Q. The actual size of the 



effect, however, turns out to be very different for the two models studied in the literature. 

An alternative approach to treat the effect of initial state radiation is by means of Monte Carlo 
event generators which include multiple gluon emission in the partem shower approximation. Compre- 
hensive phenomenological analyses have been carried out for charmonium production at the Tevatron 
and at the LHC © [@] [Q using the event generator Pythia [ ]37[ ] supplemented by the leading 



colour-octet processes [93]. The inclusion of initial state radiation as implemented in Pythia leads to 
an enhancement of the short-distance cross section. The size of the effect is significantly larger than for 



the Gaussian A; t -smearing mentioned above, and it extents out to large p t . Consequently, the (O s 



0rl 



Sn 



and (Og [ 3 -Po]) NRQCD matrix elements estimated from the Monte Carlo analysis of the Tevatron cross 
sections are significantly lower than the ones listed in Table 13 (see [93] [ f?4| ] for details)]] Figure 



shows the individual contributions to the direct J/tp cross section at the LHC as estimated with the 
Pythia Monte Carlo Note that for consistency the curves are based on the NRQCD matrix el- 



10 



10 



BR(JA)/^u. V) da(pp->JA|/+X)/dp T (nb/GeV) 

I y I < 2.5 

LHC 



color singlet 
color octet Sj 
color octet S + Pj 



10 



10 



-2 



10 



CTEQ2L 







10 



20 



30 



40 50 

P T (GeV) 



Fig. 35: Cross sections for J/ip production in pp—*.J/x}j + X at the LHC (^/~s — 14 TeV, rapidity cut \y\ < 2.5) obtained from 
a Monte Carlo event generator [p4|. CTEQ2L parton distribution functions [g^]; (i) dotted line: colour-singlet, (ii) dashed line: 
colour-octet 1 So + 3 Pj, (iii) dot-dashed line: 3 5'{ 8 ' ) , (iv) solid line: all contributions. NRQCD matrix elements as specified in 



ements extracted from the Monte Carlo analysis of the Tevatron data rather than on the values of the 

'Support is added to the Monte Carlo extraction of the NRQCD matrix elements by analyses of J/ip production in inelastic 
7p-scattering [@ || and B decays which seem to prefer small values of (Of ^So]) and {O$[ 3 P }). 



leading-order fit listed in Table |l3|. One observes that the final prediction is consistent with the result 
presented in fig. [34] within errors. The extrapolation of the Tevatron fits to LHC energies seems rather 
insensitive to the details of the underlying theoretical description, and different approaches yield similar 
predictions for the LHC cross sections as long as the appropriate NRQCD matrix elements are used. 
The Monte Carlo implementation [^] should therefore represent a convenient and reliable tool for the 
experimental simulation of quarkonium production processes at the LHC. 

A crucial test of the NRQCD approach to charmonium production at hadron colliders is the anal- 
ysis of J/ip and ip(2S) polarization at large transverse momentum. Recall that at large p t , tp production 
should be dominated by gluon fragmentation into a colour-octet 3 Si charm quark pair, fig. ]32|(c). When 
Pt S> 2m c the fragmenting gluon is effectively on-shell and transverse. The intermediate cc pair in the 
colour-octet 3 Si state inherits the gluon's transverse polarization and so does the quarkonium, because 
the emission of soft gluons during hadronization does not flip the heavy quark spin at leading order in 
the velocity expansion. Consequently, at large transverse momentum one should observe transversely 
polarized J/ip and ip(2S) 1 1 00 ] . The polarization can be measured through the angular distribution in 
the decay ip— >l + l~ , given by dT/dcos 9 oc 1 + a cos 2 9, where 9 denotes the angle between the lepton 
three-momentum in the ip rest frame and the ip three-momentum in the lab frame. Pure transverse polar- 
ization implies a = 1. Corrections to this asymptotic limit due to spin-symmetry breaking and higher 
order fragmentation contributions have been estimated to be small [101]. The dominant source of depo- 
larization comes from the colour-octet fusion diagrams, fig. [32|(d), which are important at moderate p t . 
Still, at 0(v 4 ) in the velocity expansion, the polar angle asymmetry a can be unambiguously calculated 
within NRQCD J84| ] [ 102 ] in terms of the three non-perturbative matrix elements [Table [l3|] that have 
been determined from the unpolarized cross section. In fig. 36 we display the theoretical prediction for 
a in ip(2S) production at the Tevatron as function of the ip(2S) transverse momentum. No transverse 
polarization is expected at p t ~ 5 GeV, but the angular distribution is predicted to change drastically 
as p t increases. A preliminary measurement from CDF [103] does not support this prediction, but the 
experimental errors are too large to draw definite conclusions. A similar picture emerges from the anal- 
ysis of J /if) polarization [ |104| ], where, however, the theoretical analysis is complicated by the fact that 
the data sample still includes J/ip that have not been produced directly but come from decays of higher 
excited states [105]. 

Polarization measurements are crucial to discriminate the NRQCD approach from the colour evap- 
oration model, where the cross section for a specific charmonium state is given as a universal fraction of 
the inclusive cc production cross section integrated up to the open charm threshold. In general, the as- 
sumption of a single universal long-distance factor is too restrictive. It implies a universal <t(x c )/o"( J/ip) 
ratio, which is not supported by the comparison of charmonium production in hadron-hadron and photon- 
hadron collisions. Still, since the colour evaporation model allows colour-octet charm quark pairs from 
gluon fragmentation to hadronize into charmonium, it can describe the p t distribution of the Tevatron 
data [57] [58] [ ^ ] [60]. In contrast to the NRQCD approach, however, the colour evaporation model 
predicts charmonium to be produced unpolarized. The model assumes unsuppressed gluon emission 
from the cc pair during hadronization which randomizes spin and colour. This assumption is clearly 
wrong in the heavy quark limit where spin symmetry is at work and soft gluon emission does not flip 
the heavy quark spin. Nonetheless, since the charm quark mass is not very large with respect to the 
QCD scale, the applicability of heavy quark spin symmetry to charmonium physics has to be tested by 
confronting the NRQCD polarization signature with experimental data. 

To definitely resolve the issue of quarkonium polarization, a high-statistics measurement extending 
out to large transverse momentum will be necessary. Such a measurement can be carried out at the LHC, 
where one expects a polarization pattern similar to that predicted for the Tevatron, see fig. 37. The 
absence of a substantial fraction of transverse polarization in tp production at large p t would represent 
a serious problem for the application of the NRQCD factorization approach to the charmonium system 
and might indicate that the charm quark mass is not large enough for a nonrelativistic approach to work 
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Fig. 36: Polar angle asymmetry a for ip(2S) production in pp^ip(2S)(—>fi, + fi ) + X at the Tevatron as a function of p t 



compared to preliminary data from CDF [103]. Parameter specifications as in fig. B3| NLO corrections to the fragmentation 



contribution [101] ||84] have not been included. The error band is obtained as a combination of the uncertainty (statistical only) 



in the extraction of the NRQCD matrix elements [Table [l3|] and the limiting cases that either (Of ( 1 So)) or (Of ( 3 -Po)} is set 
to zero in the linear combination extracted from the data. 



in all circumstances. 

The application of NRQCD should be on safer grounds for the bottomonium system. As v 2 ~ 0.1 
for bottomonium, higher-order terms in the velocity expansion (in particular colour-octet contributions) 
are expected to be less relevant than in the case of charmonium. Cross sections for the production of T 
states have been measured at the Tevatron in the region p t <i 20GeV J106 ] [ 107 ] [108]. The leading- 
order colour-singlet model predictions underestimate the data, the discrepancy being, however, much less 
significant than in the case of charmonium. Given the large theoretical uncertainties in the cross section 
calculation, in particular at small p t <^ My, the need for colour-octet contribution is not yet as firmly es- 
tablished as for charmonium production. The inclusion of both next-to-leading order corrections and the 
summation of soft gluon radiation is required to obtain a realistic description of the T cross section in the 
Pi -range probed by present data. Such calculations have not yet been performed, and we have therefore 
not attempted a systematic fit [ f79| ] [ p0| ] of the bottomonium NRQCD matrix elements. Our predictions 
for the T cross section at the LHC, figs. ^8 ,39, are based on a simple choice of the non-perturbative input 



parameters [Table 14 ] which is consistent with the present experimental information from the Tevatron. 
The cross sections should thus not be regarded as firm predictions of NRQCD but rather as order-of- 
magnitude estimates. The expected theoretical progress and more experimental information will allow a 
more precise prediction in the near future. 

The impact of initial state gluon radiation on the T cross sections at the Tevatron has been esti- 
mated by adding a Gaussian k t -smearing as discussed previously in the context of charmonium produc- 
tion. An average (kt) ~ 3 GeV and a K-f actor ~ 3 are found to bring the leading-order colour-singlet 




Table 14: NRQCD matrix elements for bottomonium production. The colour-singlet matrix elements are taken from the 
potential model calculation of [ |s2| [0], The colour-octet matrix elements have been determined from the CDF data for 
pt > 8 GoV Jl07| , where (Of ^So]) = (Og[ 3 P ]}/ml has been assumed for simplicity. Parameters: CTEQ5L parton 
distribution functions Mm, renormalization and factorization scale ji = {pi + Ami) 1 / 2 an d m b = 4.88 GeV. 
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cross section in line with the experimental T(15, 2S) data at p t <, My Ql09| ]. Similar results have been 
obtained within a Monte Carlo analysis [1 10], leading to significantly lower fit values for the colour-octet 
NRQCD matrix elements than those determined from a leading-order calculation [Table 0]. Moreover, 
the Monte Carlo results imply that no feeddown from x states produced through colour-octet 3 Si bb 
states is needed to describe the inclusive T cross section, in contrast to what is found at leading-order. 
The calculation of next-to-leading order corrections and a systematic treatment of soft gluon radiation 
within perturbation theory are required to resolve these issues. Figure ^ shows the inclusive T(15) 




Fig. 38: Cross sections for T(IS') production in pp—>T(lS) + X at the LHC (a/s = 14 TeV, pseudorapidity cut \r}\ < 2.5). 
Parameters: CTEQ5L parton distribution functions [|3l|], factorization and renormalization scale n — \J p'j + 4mj^, mj = 
4.88 GeV. The leading logarithms (a s lnp^ / (2m(,) 2 )™ have been summed by solving the Altarelli-Parisi evolution equations 
for the gluon fragmentation function. NRQCD matrix elements as specified in Table The error band is obtained by varying 
the colour-octet matrix elements between half and twice their central value for illustration. 



cross section at the LHC as obtained from the Monte Carlo calculation [|1 10p. The curves are based on 

As 



the NRQCD matrix elements extracted from the Monte Carlo analysis of the Tevatron data [ |11( 
in the case of charmonium production, one observes that the final LHC prediction is consistent with the 
leading-order result presented in fig. 38 within errors. 

Let us finally present the polarization pattern predicted for direct T(15) production at the LHC, 
fig. 41, based on the NRQCD matrix elements of Table 14. Higher-order corrections to the gluon frag- 
mentation function [101] [|4|] will lead to a small reduction of the transverse polarization at large p t 
and should be included once precise data become available. If the charmonium mass is indeed not large 
enough for a nonrelativistic expansion to be reliable, the onset of transverse T polarization at p t S> -Mr 
may become the single most crucial test of the NRQCD factorization approach. 

In summary, we have discussed some of the phenomenological implications of the NRQCD ap- 
proach for quarkonium production at the LHC and presented 'state-of-the-art' predictions for ip and T 
differential cross sections and polarization.^] Among the theoretical issues that need to be addressed in the 
future are the calculation of higher-order QCD corrections, the summation of higher-order terms in the 
velocity expansion and quantitative insights in the effect of higher-twist contributions. Besides a global 
analysis of different production processes and observables at various colliders, quarkonium physics at 
the LHC will play a crucial role to assess the importance of colour-octet processes and to conclusively 
test the applicability of non-relativistic QCD and heavy-quark spin symmetry to the charmonium and 



Other processes that hav e bee n studied in the literature include quarkonium pro duction in a ssociation with photons [111] 



[112] or electroweak bosons [1 13], as well as r\ and P-wave quarkonium production [1 14] [115]. 




Fig. 39: Same as fig. [?8] for T(2S) production. 



bottomonium systems. 

7. PROSPECTS FOR b PRODUCTION MEASUREMENTS AT THE LHC[] 

Of the existing and currently proposed accelerator facilities, the LHC will yield the largest rate of b 
quarks. A well defined program for b production investigations, and the development of dedicated de- 
tection strategies optimised for ATLAS, CMS and LHCb, are required for the succesful exploitation of 
the rich LHC potential. After an introduction summarising the main physics motivations, we review the 
detector and trigger features relevant for b production in the LHC experiments. The kinematic ranges 
accessible to the three experiments are then described. Theoretical motivations and possible measure- 
ment methods are presented for single b quark properties, correlations in b production, multiple heavy 
flavour production, polarization, and charge asymmetry effects in i3-hadron production in pp interac- 
tions. Based on earlier performance studies, the potential for these measurements is estimated and some 
preliminary results are presented. We conclude with a summary of the present status of the preparations 
for 6-production studies. 

7.1 Introduction 

While many LHC studies have been devoted to S-decays, b production has not yet been directly ad- 
dressed. Even though b decay investigations will provide some information on the production, at the 
discussions of this workshop it became clear that they are not sufficient to cover all aspects of produc- 
tion. 

Heavy quark production in high energy hadronic collisions is important for the study of Quantum 
Chromodynamics (QCD). Nowadays, QCD is recognized as a well established and solid theory. If 
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Fig. 40: Cross sections for T(IS') production in pp^T(lS) + X at the LHC (^fs = 14 TeV, rapidity cut \y\ < 2.5) obtained 
from a Monte Carlo event generator [110]. CTEQ2L parton distribution functions (i) dot-dashed line: 3 S{ 8 \ (ii) solid 



line: all contributions. NRQCD matrix elements as specified in [ 1 10] 



disagreements between the theoretical predictions and the experimental data are found, they will suggest 
the lack of understanding only of a particular production mechanism. In many cases these disagreements 
may be attributed to a too slow convergence of the perturbation series. In other cases, there may be 
important contributions from nonperturbative effects. Strictly speaking, the production measurements 
are not going to test the principles of QCD, but rather to outline the boundaries, where the predictions of 
perturbation theory provide an adequate description and exhaust all the visible effects. In this context, it 
will be certainly useful to test as many different processes as possible. 

We present below some examples of such processes and observables, which can potentially be 
studied in the LHC experiments. Besides testing QCD, there exist other motivations to understand pro- 
duction properties; for instance, as a control of the systematics in CP violation. Double b pair production 



is also a background in some channels of Higgs detection for LHC [116]. Measurements of the b pro- 
duction by ATLAS and CMS in the initial years of low luminosity running will also be used to optimise 
the trigger selections at high luminosity for rare B decays. 

7.2 Detector and trigger characteristics relevant for 6-production 



The ATLAS, CMS and LHCb detectors and triggers are described in detail elsewhere [ |1 17J . Even though 
the signal-to-noise ratio for b events is higher at LHC than at lower energy hadron machines, only about 
1% of the non diffractive inelastic collisions will produce 6-quark pairs. Events with B hadrons can be 
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Fig. 41: Polar angle asymmetry a for direct T(IS') production in pp— >T (IS) (— > /i + + X at the LHC as a function of pt. 
NRQCD matrix elements as specified in Table [w| other parameters as in Figure The error band reflects the limiting cases 
that either (Og^So)) or (0g'( 3 Po)} is set to zero in the linear combination extracted from the data. 



distinguished from other inelastic pp interactions by the presence of leptons, of secondary vertices and 
particles with high px- Each of the three experiments will have several levels of triggers to efficiently 
select the interesting events containing B hadrons while maintaining manageable trigger rates. The in- 
formation from the muon detectors and the electromagnetic and hadronic calorimeters will be used by the 
lowest level trigger in all the three experiments. In LHCb the lowest level trigger performs a pile-up veto 
followed by soft cuts on first level trigger objects like muons (pr > 1 GeV), electrons (Er > 2.1 GeV) 
or hadron clusters (Et > 2.4 GeV) reducing the trigger rate to 1 MHz. The more time-consuming 
operations, like vertex reconstruction and using information from the RICH for particle identification, 
will be performed by the higher level triggers. The final event rate expected from LHCb is ~ 200 Hz. 
ATLAS and CMS are central detectors for high pt physics designed to operate at high luminosities. The 
low-level trigger objects have higher pt limits than in LHCb: single muons pt > 6(7) GeV in AT- 
LAS(CMS) or dimuon triggers with a minimal pt of each muon in the interval (3 — 6) GeV in ATLAS 



and (2 — 4) GeV in CMS Ql 18| ] Ql 19Q . However, thanks to the higher luminosity, despite the higher 
Pt thresholds, they will have statistics comparable to LHCb in many exclusive channels. Simulations 
done on both experiments have demonstrated that at a luminosity of 10 33 cmr 2 s~ x , in spite of 2-3 pileup 
events on the average accompanying the b event in the same bunch crossing, i?-decays can be triggered 



on and further cleanly separated from background in off-line reconstruction [12C] [121] 



7.3 Kinematic ranges 

The central detectors ATLAS and CMS will cover the pseudorapidity region \rj\ < 2.5; the more forward 
LHCb is optimised for 1.8 < rj < 4.9 . The overlap between the experiments is less then a unit of pseu- 



dorapidity, in the region 1.8 < r] < 2.5. The low transverse momentum cutoffs in each experiment are 
limited mainly by the admissible low-level trigger rates. In the statistically dominant channels, ATLAS 
and CMS will be efficient for S-hadrons with px > 10 GeV and LHCb for px > 2 GeV. The domains 
of the Bjorken x variable for different values of the b quark transverse momentum px are given in fig. ff2| 
for two situations: when both the b and b are in a fiducial volume of a detector; and when only one of 
them is there. It is clear that in all three LHC experiments the sampled range of x is contained within 



the region already covered by HERA [122]. For comparison, the analogous distribution is calculated for 



CDF conditions (fig. 43) 
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Fig. 42: Bjorken x region of LHCb for different values of the b pr for the situation when one of the quarks is in the detector 
volume (a), both b and b are in the detector volume (b). In (c) and (d), analogous distributions are given for ATLAS/CMS. 



7.4 Single b quark production 

7.41 Theoretical motivations 

The inclusive differential cross section da/dprdrj, where pr and -q are the transverse momentum and 
the pseudorapidity of the b or b quarks, provide the basic information on b production. As discussed 
in the previous section, next-to-leading order (NLO) calculations give a cross section lower than CDF 
and DO data by a factor of ~ 2.4 [12]. However, the shape of the px distribution is well reproduced 
by LO+NLO predictions, by a semihard model of the BFKL type [|123[] and also by Pythia []124[]. In 



the region of high pr, the effects of higher order contributions are taken into account by means of the 
resummation technique [125] [23] [126]. Inref. [|23|] LO+NLO contributions are included together with 
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Fig. 43: Bjorken x region of CDF for different values of b pr for the situation when one of the quarks is in the detector volume 
(a), both b and b quarks are in the detector volume (b). 



the resummation of all terms of order a k s ln k (pT/rrib) and at ln fc (j>t / 'mj) . These contributions change 
the shape of the p? spectrum. Thus measurements of high p? single 6-spectra may be considered as a 
dedicated test for the QCD resummation technique. 



7.42 Measurement possibilities 

Experiments can measure the doubly differential cross section da/dpTdrj, where px and 7] are the trans- 
verse momentum and the pseudorapidity of a B hadron, or of a jet associated with a B hadron, or only 
one of the decay products of a B hadron (for example J/ip or //). From these experimentally measured 
quantities, the da/dpTdr] of the parent 6-quark can be extracted, using appropriate models of hadroniza- 
tion and decay. 

The determination of the absolute value of the cross section is also important. Three independent 
measurements (ATLAS, CMS and LHCb) can be done at the same energy. The determination of the 
absolute cross sections is always difficult, since it requires a precise understanding of the luminosity, 
of the trigger and reconstruction efficiency and of the background contributions. Several techniques of 



luminosity measurement are under study. It appears that precisions of ~ 3% could be achieved [127]. 
The overlap in the detection phase space of ATLAS, CMS and LHCb, in the region 1.8 < r) < 2.5 and 
Pt > 10 GeV, can be used for cross-checks. 

7.43 Exclusive channels 

From trigger and offline studies and the present experience with CDF it is known that the three LHC 
experiments can provide high statistics samples of some exclusive 5-decay channels cleanly separated 
from the background. The statistically dominant channels are those containing J/tp — > p + p~ (Bd — ► 
J/vpK , B d J/t/jK*, B ± -> J/V^ and 5° J/ip<j>), which are also needed for CP violation 
studies. Moreover, LHCb will cleanly separate large statistics of purely hadronic exclusive decays, where 
the dominant ones are Bd — > D*~ir + and Bd — * D*~af. With these processes one can cover the 
differential px cross section measurements starting approximately from pr > 10 GeV for ATLAS and 
CMS and pt > 2 GeV for LHCb respectively. The numbers of these events after three years of run at 
luminosities of 10 33 cm~ 2 s~ l for ATLAS and CMS and five years at 2 • 10 32 cm~ 2 s _1 for LHCb, are 



shown in fig. |44| as a function of a minimal transverse momentum of the 5-hadron px- 



7.44 Inclusive b — > J /ipX channels 

The inclusive channels 6 — > J/ipX can be used to extend the available statistics for production measure- 
ments to high transverse momenta (fig. 44). 



A preliminary study from CMS Q128| ] shows that for p T ~ 300 GeV, which corresponds to 
p h T ~ 550 GeV, a b-tagging efficiency of ~ 50% can be achieved with a J/ip mass and decay length 
reconstruction. This will give a signal-to-noise ratio of ~ 2.5 taking into account the prediction for 
prompt J/ip production of ref. [93]. 

In ATLAS a study has been done [ |129| ] for events 6 — * J/ip(fj,fi)X in which the pt of the b quark 
was chosen larger than 50 GeV. In particular, it was shown that the mass resolution of the J/ ip will not 
be degraded due to events in which a signal reconstructed in the muon system is wrongly associated to a 
non-muon track in the inner detector. 
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Fig. 44: Number of triggered and reconstructed B-hadrons as function of a lower cut on the B-hadron transverse momentum pr- 
Figure (a) shows the ATLAS expectations for inclusive and exclusive B-hadrons decays to J/ip after 3 years, with an integrated 
luminosity of 30ffe _1 . The full line corresponds to inclusive events, the dashed line to the sum of all exclusive channels. Fig. (b) 
shows the LHCb expectations for exclusive B-hadrons decays after 5 years. The solid line is for all statistically dominant 
channels, the dashed line shows only the channels with a J/ip in the final state. 



7.45 Inclusive b-jet production 

Another method for b production studies discussed at the workshop was based on inclusive 6-jet recon- 
struction. In both ATLAS and CMS this technique was developed for the Higgs search [13C] [ 121 ]. The 
6-jet cross section is expected to be a small fraction (close to or larger than 2% for jets with Et larger 
than about 20 GeV) of the single-jet cross section [131] [132]. If this method is to be used for single b 
quark production, it will require prescaling of the trigger for the lower px region or a cut on very high 
transverse momenta (pr > 150 GeV), to reduce the huge rate of non-6 QCD background [133]. 



Figure |45| shows the preliminary results of the CMS 6-tagging efficiency and mistagging proba- 
bility for high Et jets using the technique described in [134]. The study demonstrates that for tagging 
efficiencies of 35% — 55% the mistagging probability is better than 2% up to Et ~ 200 GeV. Beyond 
that, the 6-tagging efficiency and mistagging probability deteriorate significantly. The algorithm will be 
further optimised, possibly including lepton identification. 

The method of b cross section determination based on inclusive 6-jet identification will be heavily 
dependent on the precise understanding of the non 6-jet rejection factors. Further feasibility studies on 
this method are necessary. 
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Fig. 45: The results of CMS study of 6-tagging efficiency and mistagging probability as a function of the jet Et. The squares 
represent the results of the phase- 1 tracker, the circles those of the phase-2 tracker. 



7.5 Correlations in b production 

7.51 Theoretical motivations 

As discussed in Section ||, the overall normalisation of the production cross section, as well as the nor- 
malisation of the inclusive b spectra, remain uncertain within a factor ~ 2 because of inherent theoretical 
uncertainties. Therefore the measurement of these values does not provide a sringent test of NLO con- 
tributions. The expected correlations between the b and the b quarks can be computed in leading and 
next-to-leading order The shapes of two-particle distributions are sensitive to the NLO contribution, 
and thus can be used for these tests. In particular, distributions in the following quantities involving both 
the b and b quark can be considered: the relative azimuthal distance A<p(bb) < 1, the pair invariant mass, 
the pair transverse momentum and the pair rapidity. 

7.52 Measurement possibilities 

The choice of the decay channels is driven by the requirement that the acceptance should not vanish when 
the b and the b are close in phase space. The goal is to avoid isolation cuts in both trigger and offline 
algorithms requiring a large separation between the decay products of the B and of the B. The processes 
under consideration are based on the reconstruction of a J/ip originating from the displaced vertex of a 
S-hadron, and of an additional lepton coming from the semileptonic decay of the associated B hadron. 
For example, in the ATLAS experiment, for an integrated luminosity of 30 ffe -1 , approximately ~ 5 • 10 5 



such events are expected, with the exclusively reconstructed S-decays containing the J/ip (Table |15|) 



CDF and DO measurements showed that bb pairs are mostly produced back-to-back [135]. However 



the region most sensitive to differences between the models is A(f>(bb) < 1 rad, where only ~ 14% of 



the events are expected [124]. The statistics may possibly be increased using the semi-inclusive decays 
66 — > J/ipX accompanied by a lepton (Table 15). As an example we quote recent studies in ATLAS 
[129], performed using simulated events with — > J/t/jK°, J/ip — > /i + . They indicate that the 
signal events can indeed be reconstructed in cases when the difference of azimuthal angles between the 



J/ip and the other muon is small (Fig. 46). It is important to note that no selection cuts requiring model 
dependent corrections were necessary. 

The study can be extended to events with J/ip —> /.i + ^~accompanied by an electron and for 
J/ip — > e + e~ combined with a muon or an electron. Using all these combinations of leptons will 
allow the measurement of the same variables by different detectors, leading to an improved control of 
systematic errors. 
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Fig. 46: Reconstruction of bb — ► J /ip([i[i)X + fi using combined Muon-Inner detector off-line reconstruction in ATLAS. 



Table 15: Semi-inclusive and exclusive channels candidates for bb studies. Statistics are given for ATLAS after 3 years, with 
an integrated luminosity of 30 f6 _1 . 



Inclusive channels 


Statistics 


Exclusive channels 
with the same 
lepton content 


Statistics 


bb -> J/tl)(nn)X + fi 
bb -> J/il)(nn)X + e 
bb -> J/i{j(ee)X + \i 


2.8 • 10 y 
3.6 • 10 6 
0.6 • 10 6 


bb — > hadJ/ifj — * /-i + [J,~ + // 
bb — * hadJ/ip — » H> + + e 
66 — > hadJ /ip — > e + e~ + // 


2.1 • 10 5 
2.1 • 10 5 
0.9 • 10 5 



7.6 Multiple heavy quark pair production 

7.61 Theoretical motivations 

At present, only the leading order calculation, O(af), is available [136] for the bbbb production cross 



section. The effects of higher order corrections can only be estimated using the event generator Pythia 
5.7. Since the predictions of PYTHIA appear to be about a factor of 10 above the leading order analytical 
calculations [ 137| ], further theoretical studies are needed. 



7.62 Measurement possibilities 

Events with four b quarks can be identified in several ways, the most appropriate one depending upon 
the context. As a background to Higgs search, the requirement is four b jets in the fiducial volume. For 
the purpose of testing QCD predictions on double b production, it may be sufficient to reconstruct events 
with three b quarks in the fiducial volume. For three b quarks wither > 10 GeV and \rj\ < 2.5, Pythia 
gives a cross section of 140 nb, which corresponds to 140 events produced per second. Despite this large 
number, it will be necessary to define features allowing on-line selection of these events in the presence 
of huge non-6 and single b backgrounds. 

As a source of an incorrect tag in CP violation measurements, the relevant bbbb events are those 
with two b hadrons, produced with the same flavour charge, identified in the fiducial volume, while two 
other 73-hadrons, produced with the opposite flavour, are not detected. A direct measurement of this 
case could use reconstructed charged mesons or baryons, which are self-tagging. However the expected 
statistics of these events is insufficient. In fact, double b production is expected to be only a minor source 



of wrong tags Q138Q . Techniques exist to determine the wrong-tag rate from all processes regardless of 



its origin, which does not need to be identified [ 138| ] 



Similar to the case of double b production, the production of doubly heavy hadrons, such as the 
B c + b + c, Eb c + b + c, etc., refers to an O(af) lowest order QCD process, and also provides a test 
of perturbative QCD calculations. The question of higher order QCD contributions and, probably, non- 



perturbative contributions is still open Q139| ]. The total production rate in this case will not be indicative 
enough to establish the role of different production mechanisms, and the measurement should instead 
concentrate on the specific event topologies. In particular, various correlations between particles carry- 
ing charm and bottom may be of importance. 

Measurement possibilities are under investigation for the channel — > J/ipir [ ]140| ] . A list 



of possible semileptonic and nonleptonic B c decays can be found for instance in Q141| ]. The decays 
B ( * ] -> J/^pfiu, B { * ] -> J/V>P + , B ( * ] -> J/ipK + *, B { * ] -» J/ipD+ are other potentially interesting 
modes. 



7.7 Other measurements 

B hadrons with non-zero spin can be polarized perpendicular to their production plane. Polarization 
measurements of b hadrons produced in nucleon fragmentation could clarify the problems of different 



polarization models [|142J] that failed to reproduce the existing data on strange hyperon production [ |143[ ]. 
In particular, information about the quark mass dependence of polarization effects could be obtained. 
For symmetry reasons, in pp collisions this polarization vanishes at zero rapidity, so that the expected 
observed polarization in ATLAS and CMS will be smaller than in the more forward LHCb. Using the 
method of helicity analyses of cascade decay A° — > A°J/ip the A° polarization can be measured in 



ATLAS with a precision better than 0.016 [144]. Another approach to A? polarisation measurement, 



using the same decay channel can be found in [J145J]. 

In proton-proton collisions a charge production asymmetry of b hadrons is expected. The asym- 
metry is defined as the difference of production probabilities of a B hadron and its antiparticle. From 
the theoretical point of view, the asymmetries can provide information on the effects of soft dynamics 



during the fragmentation and hadronization (i.e., on the soft interactions between the produced b quark 
and the remnants of the disrupted proton). The relevant physical effects are expected to be unimportant 
[ |47| ] [146] [[M]] in the central rapidity region covered by ATLAS and CMS. In the more forward region 
of LHCb the asymmetry may rise to a few percent. A detailed theoretical discussion of this issues are 
given in Section |5[ 

Any production asymmetry is always measured in the presence of a CP violation asymmetry orig- 
inating from S-hadron decays. In some cases these two effects are expected to be of the same order. This 
is the case, for instance, in the channels Bd — > J/ipK*(K + Tr~), B + — > J/ipK + and A° — ► A°J/xp, 
which are expected to have a small CP violation (< 1%). A way of estimating the relative size of these 
two effects may be based on the fact that the production asymmetry varies with the transverse momentum 
and the rapidity of produced 6-quark, while the decay asymmetry should remain the same. Measurements 
of such small effects will require good understanding of the possible instrumental detection asymmetries. 



7.8 Conclusions 

The properties of b production at the LHC can be measured by the three experiments, which are com- 
plementary in phase space. The small overlap region will allow a cross check on the cross section 
normalization. The kinematic conditions are such that Bjorken x values sampled in b-production are 
above 10~ 5 , a region lower than at the Tevatron, but already covered by HERA. Differential cross sec- 
tion measurements using exclusive B hadron decays will be most important at small pt values. At high 
Pt values and for correlations and multiple heavy flavour production measurements, the statistics can 
be increased by semi-inclusive B decays containing J /tp. Possible methods using 6-jets require further 
study, to control the non-6 QCD background. The enormous LHC statistics will also allow to study the 
production polarization and charge production asymmetries. 



8. TUNING OF MULTIPLE INTERACTIONS GENERATED BY PYTHIAQ 
8.1 Introduction 

The track multiplicity distribution as well as the transverse momentum distribution of charged particles in 
proton-proton interactions (the so-called minimum bias events) affect the performances of the low level 
triggers and the detector occupancy of the LHCb experiment [147]- They should therefore be modelled 
reliably in Monte Carlo programs. In particular, at LHC energies, multiple interactions play an important 
role, and should not be neglected. 

In Section $2 we examine the multiple interaction models available in Pythia [|7|] to describe 
the event structure in hadron-hadron collisions. In Section 8.3 we select a compilation of homogeneous 
data at different energies suitable to tune the multiple interactions parameters of Pythia; the tuning 
procedure is presented in Section 8.4. We use the phenomenological extrapolations at LHC energy 
in order to get the predictions for the track multiplicity and the transverse momentum distributions in 



minimum bias and bb events; these are reported in Section 8.5 



8.2 Multiple interaction models 

The multiple interactions scenario is needed to describe the multiplicity observables at hadron collid- 
ers [148] and is also supported by direct observation [149] [ ]150| ] Jl 51 ]- The basic assumption is that 
several parton-parton interactions can occur within a single hadron-hadron collision. Four different mod- 
els are available in Pythia. The main parameter of these models, Pt , is the minimum transverse 
momentum of the parton-parton collisions; it effectively controls the average number of parton-parton 
interactions and hence the average track multiplicity. The differences between the four models, which 
mainly affect the shape of the multiplicity distribution, are the following: 
10 Section coordinators: P. Bartalini, O. Schneider. 



Model 1 (default in Pythia) 

All the hadron collisions are equivalent (as opposed to model 3 and 4 below) and all the 
parton-parton interactions are independent; the Pr min parameter represents an abrupt cut-off. 
Model 2 

Same as Model 1 but with a continuous turn-off of the cross section at Pr ■ . 
Model 3 

Same as Model 2, but hadronic matter in the colliding hadrons is distributed according to a 
Gaussian shape, and a varying impact parameter between the two hadrons is assumed. 
Model 4 

Same as Model 3 but the matter distribution is described by two concentric Gaussian distribu- 
tions. 

The varying impact parameter models (Models 3 and 4) were developed []148|] to fit the UA5 



data [152]. A recent study performed by the CDF collaboration [153] concludes that a varying impact 
parameter model (Model 3) is also preferred to describe the underlying tracks in b events produced at the 
Tevatron. 

In the absence of published results on multiplicity distributions in minimum bias events at the 
Tevatron, we compare again the predictions of Model 1 and Model 3 with the UA5 data, using the final 



charged multiplicity distribution from the full pp data sample collected by UA5 at \fs = 546 GeV [154], 



a recent version of PYTHIA^], and a modern set of parton distribution functions, CTEQ4L [14], tuned on 
both HERA and Tevatron data. For this comparison, the main multiple interactions parameter (Pr min ) 
is tuned in each model to reproduce the mean value of the measured charged multiplicity distribution in 
not single diffractive eventsQ (< N ch >= 29 A ± 0.3 ± 0.9). We obtain P Tmin = 1.63 ± 0.02 for Model 
1 and Pr min = 1-97 ± 0.03 for Model 3. The shapes of the multiplicity distributions are compared in 
fig. It is clear that Model 3 is preferred over Model 1 to describe the UA5 data. In particular the shape 
of the tail at high multiplicities is reproduced much better by Model 3. The UA5 results are corrected for 
the lower efficiency expected on double diffractive events. Therefore the simulated samples include the 
generation of all kind of non single-diffractive events. The uncertainty in the diffractive cross sections 
relative to the partonic ones can affect the observed discrepancies between the data and the PYTHIA 
predictions in the low multiplicity region^. 

8.3 Mean charged multiplicity at t] = 

In order to produce realistic Pythia predictions for the multiplicity observables in the LHC environment, 
it is necessary to take into account the energy dependence of the Pr min parameter. Unfortunately there 
are not many published data concerning the charged multiplicity distribution in minimum bias events at 
hadron colliders. On the other hand there are some data available relative to the average charged mul- 
tiplicity in non single-diffractive events, in particular for the central pseudo-rapidity region. Therefore, 
to study the energy dependence of the Pr min parameter at generator level, we consider an homogeneous 
sample of corrected average charged multiplicity measurements at six different center-of-mass energies 



(y^ = 50, 200, 546, 630, 900 and 1800 GeV) in the pseudo-rapidity region |r/| < 0.25 [156] [[1571 



The energy dependence of dN c h/drj at rj = is shown in fig. |48| a together with the fit of a quadratic 



function of ln(s) proposed in Reference [157]; using this fit to extrapolate at LHC energy would predict 
dN ch /drj = 6.11 ± 0.29 at rj = 0. 



"Version 6.125 [|155j] was used for all the studies reported here. 

12 In this paper we define as "non single-diffractive event" any inelastic hadron-hadron interaction that cannot be regarded as 
a single diffractive event; in the framework of the PYTHIA hadronic interactions, the "non single-diffractive" sample includes 
the 2 — + 2 partonic processes and the double diffractive hadron-hadron interactions. 

13 The pp cross sections predicted by PYTHIA at y/s = 546 GeV are 30.7 mb for the partonic processes and 5.3 mb for the 
double diffractive processes. 
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Fig. 47: Charged multiplicity distribution for non single-diffractive events in pp collisions at 



546 GeV as measured 



by UA5 [154] compared with PYTHIA predictions using the CTEQ4L parton distribution functions and either Model 1 (solid) 



or 3 (dashed) for multiple interactions. In each case the Pr min parameter has been tuned to reproduce the mean multiplicity 
measured in the data. 



8.4 Tuning of the multiple interaction parameter Pt 

The average charged multiplicity measurements performed on non single-diffractive data in pp collisions 
and described in Section |83| are used to tune the main multiple interaction parameter in PYTHIA, Pr min ■ 
We generate non single-diffractive events. At each value of yfs, the Pr min parameter is adjusted to 
reproduce the average multiplicity measured in the data. The uncertainty on the tuned value of Pr min 
reflects the uncertainty on the data. However, the tuned parameters depend on other aspects of the 
PYTHIA simulation: in particular the effects of various choices for the multiple interaction model and 
the parton distribution functions are investigated. For simplicity, the results of these studies are shown 
only for some representative settings: 



as an example of pre-HERA parton distribution functions we consider the CTEQ2L Q99| ] set used 
by default in PYTHIA versions 5.7, but recently retracted by their authors; 



as an example of post-HERA parton distribution functions we consider the CTEQ4L [14] and 



GRV94L [ 158 ] sets, the latter being the new default in PYTHIA versions 6. 1. 



This study is restricted to Models 1 and 3 for multiple interactions (see Section 8.2). 

The results of the tuning procedure are shown in fig. |48fc >: in each case, Pr min appears to be 
monotonically increasing as a function of y/s. This dependence is much more pronounced for the post- 
HERA parton distribution functions regardless of the choice for the multiple interactions model. 

It was shown in Reference fll59j ] that the post-HERA parton distribution functions imply an 
energy-dependent Pt cut-off. This is heuristically motivated by the Lipatov-like dependence of the 
gluonic parton distribution function in the small-x limit: 

xg(x, Q 2 ) — ► constant x x~ € for x — > (35) 

with e ~ 0.08, while the pre-HERA parton distribution functions give a reduced charge screening effect 
and consequently a less sensitive running of the Pt cut-off. This is heuristically motivated by the Regge- 
like dependence of the gluonic parton distribution function in the small-x limit: 



xg(x,Q 2 ) — > constant for x — > 



(36) 
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Fig. 48: a) dN c h /dtq at r\ = as a function of y/s. The solid curve represent a phenomenological fit performed by CDF [ 157 ] 
with the formula dN ch /dr](s) = (0.023 ± 0.008) ln 2 (s) - (0.25 ±0.19) ln(s) + (2.5 ± 1.0). The two dashed curves represent 
the 1 sigma variations of the fitted parameters. 

b) y^s-dependence of the Pr min parameter for various parton distribution functions and multiple interactions models. The 
points with error bars are the results of the tuning procedure on the data. The curves are the results of the fits through the points 
assuming the functional form of Equation ^7] and are characterized by the parameters given in Table [Ifj. 



Table 16: Results of the fits describing the exponential running of the PYTHIA multiple interactions parameter Pr min for 
different parton distribution functions and multiple interactions models. 



Mult. int. model 


PDF 


P^ HtJ (GeV/c) 

rn.m. 


e 


3 


CTEQ2L 


1.99 ±0.11 


0.048 ± 0.006 


3 


GRV94L 


4.06 ± 0.24 


0.103 ± 0.006 


3 


CTEQ4L 


3.47 ±0.17 


0.087 ±0.005 


1 


CTEQ4L 


3.12 ±0.16 


0.100 ± 0.005 



In order to extrapolate Pr min at LHC energy, one needs to find a reasonable function to fit the tuned 
Pr min values as a function of y/s for the different parton distribution functions and multiple interactions 
models; a four degree of freedom fit is performed using the following exponential form, inspired by the 



recent implementations added in PYTHIA since version 6.120 [ |155| ]: 



PT min {yTs) = P : 



LHC 



14TeV 



2c 



(37) 



The fitted functions are superimposed on fig. 48b and the results obtained for the fitted parameters e 



and HC are given in Table [l6[ This quantitative analysis demonstrates that the power law expressed in 
Equation 37 holds for values of e between ~ 0.08 and ~ 0.10 if post-HERA parton distribution functions 
are used, and for somewhat smaller values of e (~ 0.05) for the pre-HERA parton distribution functions. 



8.5 Pythia predictions at LHC energy 



pLHC 

-^777 i n. 



PiTevatron 
T ■ 

J mm 



2c 



(38) 



Figures |49| a and b show multiplicity and pseudorapidity distributions for charged particles predicted by 
Pythia at LHC with CTEQ4L and Model 3 for multiple interactions. The value of -Pr m „ at LHC is 
obtained by an extrapolation 

14 TeV 
1.8 TeV 

where pTevatron j s ^ e p value tuned at the Tevatron energy of 1.8 TeV. For the parameter e, we 
adopt the results given in Table [lg. It is important to note that the predictions dN^/dr} = 6.30 ± 0.42 
(for e = 0.087 =F 0.005) at rj 
(dN ch /dr] = 6.11 ±0.29). 

In order to demonstrate the importance of the correct Pp, 



are consistent with the phenomenological fit displayed in fig. 48 a 



results obtained by assuming P^ HC 



extrapolation, figs. |49| a and b also show 
i.e. e = not supported by the data as demonstrated 
in Section The multiplicity distribution has a tail at high multiplicities and dN c h /drj at r\ = is not 
consistent with that obtained from the phenomenological fit. 
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Fig. 49: PYTHIA predictions at LHC energy, using the CTEQ4L parton distribution functions and Model 3 for multiple 
interactions with Pr min given by Equation the solid and dashed distributions correspond to the central value and ±1 a 
uncertainties of the fitted value e = 0.087 ± 0.005. The dotted histogram is obtained with e = 0, i.e. using the Pr min value 
tuned on Tevatron data and ignoring its energy dependence, a) Charged track multiplicity in the entire 4n solid angle, b) 
Average charged track multiplicity per unit pseudorapidity, as a function of pseudorapidity. 



Figure 50 shows the same distributions as fig. but for the CTEQ2L parton distribution func- 
tions. It is interesting to note that, once the extrapolation of Pr min is properly done, there is no large 
difference between the multiplicity and pseudorapidity distributions obtained with different structure 
functions. 

Figure [5T|a-d compare the Pythia predictions for the multiplicity and transverse momentum dis- 
tributions in the LHCb angular acceptance (1.8 < rj < 4.9) between minimum bias and bb events^. 
These predictions are obtained with CTEQ4L, multiple interactions Model 3 and the proper Pr min ex- 
trapolation. They show clear differences between minimum bias and bb events, in particular higher 
average multiplicity and transverse momentum for bb events. 

14 The bb events are selected among the minimum bias events. 
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Fig. 50: PYTHIA predictions at LHC energy, using the CTEQ2L parton distribution functions and Model 3 for multiple 
interactions with Pr min given by Equation the solid and dashed distributions correspond to the central value and ±1 a 
uncertainties of the fitted value e = 0.048 ± 0.006. The dotted histogram is obtained with e = 0, i.e. using the Pr min value 
tuned on Tevatron data and ignoring its energy dependence, a) Charged track multiplicity in the entire 4-zr solid angle, b) 
Average charged track multiplicity per unit pseudorapidity, as a function of pseudorapidity. 



Pythia predictions with multiple interactions Model 1 for the multiplicity and transverse momen- 
tum distributions are shown in fig. 52 for minimum bias and bb events. Compared to the results obtained 
with Model 3, less significant differences between minimum bias and bb events is observed. In Sec- 
tion |8]2| we have stressed that a varying impact parameter model for multiple interactions (i.e. Model 3) 
is needed to describe the charged track multiplicity in hadron-hadron interactions. There are arguments 
in favour of adopting a multiple interactions model with varying impact parameter to describe the heavy 
flavour production at hadron colliders [160], though there are no experimental data at low transverse 
momentum and high pseudorapidity (i.e. in the LHCb acceptance region). A more detailed discussion 
on the effect of multiparton interactions in bb events at LHCb can be found in Reference [161]. 



8.6 Conclusions 

Comparisons between Pythia and experimental data demonstrate that, in order to reproduce the charged 
track multiplicity spectrum, a varying impact parameter model has to be adopted. 

The varying impact parameter models predict sensitive differences in multiplicity and Pt distri- 
bution between light and heavy flavour events. 

The running of the Pt cut-off parameter in PYTHIA multiple interactions is mandatory. Predic- 
tions made at LHC energy with a fixed Pt cut-off tuned at lower energies overestimate the multiplicity 
observables. Taking into account the running of the Pt parameter is even more important if post-HERA 
parton distribution functions are used. 
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Fig. 51: PYTHIA predictions for charged tracks in the LHCb acceptance using the CTEQ4L parton distribution functions and 
Model 3 for multiple interactions with proper Pr min ', the normalized predictions for bb events (dashed curve) and minimum 
bias events (solid curve) are superimposed, a) Charged track multiplicity distribution, b) Transverse momentum distribution of 
charged tracks, c) Distribution of the maximum transverse momentum of charged tracks, d) Rejection efficiency as a function 
of a traverse momentum cut (an event is rejected if all the charged tracks have a transverse momentum below the cut). 
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